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Abstract 

We study the critical two-dimensional Ising model with a defect 
line (altered bond strength along a line) in the continuum limit. By 
folding the system at the defect line, the problem is mapped to a 
special case of the critical Ashkin- Teller model, the continuum limit 
of which is the Z2 orbifold of the free boson, with a boundary. Possible 
boundary states on the Z2 orbifold theory are explored, and a special 
case is applied to the Ising defect problem. We find the complete 
spectrum of boundary operators, exact two-point correlation functions 
and the universal term in the free energy of the defect line for arbitrary 
strength of the defect. We also find a new universality class of defect 
lines. It is conjectured that we have found all the possible universality 
classes of defect lines in the Ising model. Relative stabilities among 
the defect universality classes are discussed. 
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1 Introduction 



Boundary (surface) critical phenomena in two-dimensional systems have 
attracted a lot of interest. Even if the bulk system is a free system, the pres- 
ence of the boundary makes the problem non-trivial. Boundary conformal 
field theory |@] is a powerful tool to solve these problems. Its applications to 
quantum impurity problems are also discussed 0. A variation of the prob- 
lem has a defect line inside a two-dimensional system. This kind of problem 
can still be mapped to a system with boundary, by folding the system at the 
defect line [|[. 

The two-dimensional Ising model, which was solved 50 years ago by On- 
sager 0, is still useful to test new techniques and ideas in statistical me- 
chanics. In this paper, we apply boundary conformal field theory to the 
two-dimensional Ising model with a defect line. This problem has been stud- 
ied by several authors. Bariev p found a continuously varying "surface" 
critical exponent. McCoy and Perk ^ and Kadanoff 0] discussed the corre- 
lation functions along the defect line. Brown discussed several properties 
of this model. In particular, he calculated two-point spin correlation function 
for general locations, in first-order perturbation of the defect strength. Sev- 
eral other studies are focused on the finite-size scaling of the transfer-matrix 
spectrum on a cylinder. Cabrera and Julien [|I4| calculated the spectrum nu- 



merically. Following the application of conformal invariance by Turban []10 



see also Ref. [inj), Henkel et. al. [T^, |T3[ studied the spectrum of the quan- 



tum version by various methods and discussed the spectrum in terms of 
Virasoro and Kac-Moody algebras. Exact formulae for the spectrum on the 
square lattice with general anisotropy are given by Abraham et. al. [p^ 



Grimm |T6[ also examined generalizations of the defect line. Recently, an 
S'-matrix approach to the problem was also presented by Delfino et. al. [|17 



In this paper, we study the continuum limit at the critical point. To relate 
this problem to boundary conformal field theory, we "fold" the system into 
a c = 1 conformal field theory with a boundary. Using the exact partition 
function, we identify the boundary states in terms of c = 1 conformal field 
theory. Folding the Ising model gives two decoupled Ising models in the bulk, 
which is a special case of the Ashkin- Teller model. The critical theory of the 
Ashkin- Teller model is a c = 1 conformal field theory: the Z2 orbifold of a 
free boson. We explore the possible boundary states of the orbifold theory 
and apply a special case to the Ising defect problem. 
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Identification of the boundary states enables us to give the complete spec- 
trum of "surface" critical exponents for general strength of the defect. Cardy 
and Lewellen |T^ showed that the boundary correlation functions can be 



determined from the boundary state and the Operator Product Expansion 
(OPE) in the bulk. We employ this idea to calculate the exact two-point spin 
correlation function, using Zamolodchikov's solution for c = 1 conformal 
blocks of spin operators. We also calculate certain four-point correlation 
functions, as well as some correlation function near the end of the defect 
line. A universal term in the free energy with a finite length defect is also 
given. Moreover, we find a new one-parameter universality class of defect 
lines in terms of boundary conformal field theory. We calculate two-point 
spin correlation functions and then identify a corresponding defect in the 
quantum lattice model. Finally, the relative stability of the defect universal- 
ity classes is discussed, and a conjecture is made regarding the complete set 
of defect universality classes. 

The organization of the paper is as follows. In Section ^, we derive the 
partition function of the model in the critical continuum limit. We use the 
result to identify the boundary condition in terms of Ashkin- Teller model in 
Section ^. The possible boundary states on the Z2 orbifold theory is explored 
in Section ^, and further applied to the present problem. In Section |], we 
calculate two-point correlation functions and other universal quantities for 
arbitrary strength of the defect. A new one-parameter family of defect lines 
is studied in Section In Section ^ we discuss the relative stability and 
renormalization-group fiow among the defect universality classes. We give a 
summary and some discussion in Section 



A brief description of the present work has appeared in [18 



2 Partition function at the critical point 

Here we derive the partition function of the Ising model with a defect line, 
in the continuum limit. While Henkel et. al. |T^, |T^ have found the result 
in the r-continuum limit (quantum Ising chain), here we derive the partition 
function for the square lattice Ising model with an arbitrary anisotropy, tak- 
ing the continuum limit of the analysis by Abraham et. al. [^. Apart from 
the trivial non-universal parts, the result is independent of the anisotropy, 
under an appropriate rescaling. It agrees with that of Henkel et. al. which 
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corresponds to the extremely anisotropic limit. This is a manifestation of 
the universality in the Ising model in the presence of a defect line. 

We consider a square-lattice Ising model with the lattice constant a. The 
Ising model on a cylinder with a defect line is defined by the classical Hamil- 
tonian: 

i=l,M-l j=l,N 3=l,N 

(2.1) 

The altered link (defect) is placed between i = M and i = 1. We denote 
the defect strength hy b = J/J2. This model reduces to the periodic bound- 
ary condition, the free boundary condition and the antiperiodic boundary 
condition, respectively for 6 = 1, 0, —1. We also introduce Ki, K2 and by 

Ki = Ji/kT, K2 = Ji/kT (2.2) 

sinhir*sinhJs:i = 1. (2.3) 

Let us refer to the direction parallel to the defect line as horizontal. 
For periodic or antiperiodic boundary condition the transfer matrix can be 
mapped to a free-fermion Hamiltonian |^ by the Jordan- Wigner transfor- 



mation. However, it is not straightforward for general b because the system 
lacks translation invariance in the i (vertical) direction and a simple Fourier 



transformation is not useful. Nevertheless, Abraham, Ko and Svrakic []15 
obtained the spectrum of the transfer matrix in the j (horizontal) direction 
using a spinor approach. We start from their result and apply it to the 
critical case. 

According to their result, the transfer matrix V can be simply expressed 
by fermion operators, but the Hilbert space is divided into two sectors: 

V = P+V+ + P-V-. (2.4) 

Here P+ {P-) is the projection operator onto the even (odd) fermion-number 
sector. However, one must be careful about the definition of the "even" and 
"odd" sectors. (See below.) The transfer matrices in the two sectors are 
given by 

M 

V± = exp[-^7K^)(c]c, - 1/2)]. (2.5) 
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The "single-fermion energy" 7 is given by the Onsager dispersion function 



cosh7(p) = cosh2i^i* cosh 2/^2 — sinh2Ki* sinh2ii'2 cosp (2-6) 
and the quantization condition of the wave number Q < k < -k / a 

tan^ ^ ^inh[(l-a:)ir2] ^p) (2.8) 
2 sinh [(1 + a;)is:2] 2 ^ ^ 

where ± in the first equation corresponds to V±, a is the lattice constant, 
( = +1 or —1 which is independent of V±. 6* is defined as 

sinh 7(p) cos S*{p) = sinh 2K2 cosh 2KI — cosh 2K2 sinh 2KI cos p. (2.9) 

Let us discuss the critical point where K2 = but b still remains as a 
free parameter. In order to discuss the universal behavior at large distances, 
we take the continuum limit a —>■ (or equivalently consider the scaling limit 
where the length scale is much larger than the lattice constant), keeping the 
circumference of the cylinder L = Ma constant. 

Then only the low-energy limit of the dispersion relation is relevant in 
the discussion. It is given by 

^{k) = v\kl (2.10) 

where f is a constant (spin wave velocity). In the following, we rescale the 
horizontal ("imaginary time") direction so that v = 1. The transfer matrix 
can be described by a Hamiltonian in the continuum limit. The Hamiltonians 
for fermion number even/odd sector are given by 

Hi = -1/2). (2.11) 

j 

In this limit, S*{k) = 7r/2 independent of k and the quantization con- 
dition ( |2.8|) is much simplified. Furthermore, we can choose freely the sign 
of k in the discussion of the energy, because the dispersion function is an 
even function of k. Thus the quantization condition in the continuum limit 
is written as 

,_L «(=F&) 27r 
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where n is an arbitrary integer and a is now independent of k: 

^''''^ - smh[{l + x)K,y ^^-^^^ 

From this equation, we see that a{x) satisfies 

a{x) + a{-x) = n. (2.14) 

a{x) takes the values vr, tt/2 and respectively for x = —1,0 and 1. When 
X ±oo, a is given by 

tan^i|^ = -e^^^^ (2.15) 

Thus 71 < a(— oo) < 3n/2 and —n/2 < a(oo) < 0. We note that, in the 
K2 —>■ limit with the x —>■ ±00 limit, a approaches the limiting values 
a(— 00) = 37r/2 and a(oo) = — 7r/2. The K2 — > limit corresponds to the 
anisotropic limit where the horizontal link becomes very strong. (It should 
be remembered that we stay at the critical point where K2 = K^.) Regarding 
the horizontal direction as imaginary time, this is the so-called r-continuum 



limit |22| of the Ising model. In this limit, the system is equivalent to the 
one-dimensional quantum Ising model in a transverse field. The quantum 
Hamiltonian is given by 

00 

H = - ^""(n) -Y.a'{n-l)a'{n) -ba'{-l)a'{0), (2.16) 

n=— 00 n^O 

where a^'^{n) is a Pauli spin operator at site n. This model is the critical 
transverse Ising model with a defective link between n = —1 and 0. 

In eq. ( |2.11| ), we arrived at a rather simple description of the transfer 
matrix for generalized boundary condition with parameter b. Namely, it 
is always a free fermion Hamiltonian and the dispersion relation ( |2.1(j| ) is 
independent of b. The only effect of the boundary condition is the "phase 
shift" a of the wavenumber as in (|2.12|). This phase shift is different in even- 



and odd- fermion number sectors. 

We note that the parity of the fermion number is a subtle problem. 
Exchanging a fermion annihilation operator Ck and the corresponding cre- 
ation operator cj, preserves all the fermion anticommutation relations. How- 
ever, this procedure flips the fermion number parity operator deflned by 
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(— 1)'^ = (— Thus the fermion parity operator depends on the def- 
inition of the fermion creation/annihilation operators. We fix the fermion 
operators so that all particles have non-negative energy and the Hamiltonian 
takes the form of eq. ( |2.11| ). For a > 0, the two sectors of Hilbert space 
in ( p.4|) actually corresponds to the even/odd fermion number in the above 
definition: 

P+ = \ ^ , P~ = (2.17) 

When a changes sign, one of the quantized wavevectors passes through k = 0. 
The energy of the corresponding fermion seems to change non-analytically 
as \k\. However, this should be understood as the fermion energy depending 
on on the parameter linearly but the annihilation/creation operators being 
exchanged when the single fermion energy becomes negative. Hence the 
fermion number parity should be flipped when a becomes negative (as long 
as we keep the above definition of the fermion operators.) 

Let us consider the partition function of the Ising model on a cylinder 
with circumference L and length (3 (a macroscopic length scale), as shown in 
Fig. |I]. We apply the periodic boundary condition in the f3 direction. The 
partition function is given by 

Z = TTexp{-f3H) (2.18) 

where H = + P H is the total Hamiltonian. 

First we consider the — 1 < & < 1 case. Here a{±b) > and thus the 
assignment ( p.l7|) is valid. In the following we denote a{b) by simply a: 



< a < 71 for —l<b<l. In the odd ((—1)'^ = —1) sector, the fermion 
one-particle energies are given by e„ = |(a + 27m) /L\. The ground state is 
the fermionic vacuum |vac) which satisfies Cfe|vac) = 0. The vacuum still has 
"zero-point energy" 

^vac = -^Een (2.19) 



according to eq. (|2.11|) . While this quantity apparently diverges (but of 



course is cut off by the lattice) in the continuum limit, we can extract the 
universal "Casimir energy" in the odd sector as 

' - ?1 (2.20) 
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where C is a non-universal energy density, which is the order of inverse 
square of the lattice constant. The second term is a non- universal shift of 
the ground-state energy due to the defect, ("surface energy"). The last term 
is the universal Casimir energy. 

For even sector, we can simply replace a by vr — a to obtain 



K1 



CL + A + 



L 



TT 

12 



a 



(2.21) 



Since there must be at least one fermion in the odd sector, the ground state 
energy in the odd sector is given by 



TT — a 



(2.22) 



The ground state energy in the even sector is simply (|2.21| ) and is lower than 
E~ . Thus eq. (|2.21|) also gives the ground-state energy of the present system. 
It reduces to the known value for periodic (a = 0), antiperiodic (a = tt) and 
free (a = 7r/2) cases [29]. 

In the following we ignore the non-universal part and set L = 1 for 
simplicity. (The L dependence can be recovered by simply replacing (3 — > 
/?/L.) The partition function in the odd sector is given by 



?odd = e"^^^^^ P^e-^T.u'^^*' (2.23) 

{nfc=0,l} 



Using 



l-(-l)^ 1_(_1)E."'= 



--2 2 ' (2-2^) 

we obtain the infinite product representation of the elliptic theta functions 
as 



^odd 



ri{q) 



ri{q) 



(2.25) 



where the parameters is given by 



q = e-'^^^,w 



-a(3 



(2.26) 
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In this paper we use the elhptic theta functions and Dedekind eta function 



di{w,q) 
ri(q) 



2*." 2> UJ^ 2 



n=—oo 
oo 

E ^ 

n=— oo 
oo 

71= — oo 
OO 



J2 i-l^q-^^' w' 



(2.27) 
(2.28) 
(2.29) 
(2.30) 
(2.31) 



n=l 



Theta functions with a single argument ^j{q) should be understood as '^j{w 
l,q) ("theta constants"). 

For the even sector we make a similar calculation and obtain 



g-a2/3/4^ 



^3{w,q) _^ ^4{w,q) 



r]{q) 



Thus the total partition function Z = Zeven + -^odd is given by 



g-a2^/47r 



^Ising 



^3{w, q) i^4{w,q) i^2{w,q) .i^i{w, q) 



(2.32) 



(2.33) 



r]{q) r]{q) r]{q) r]{q) 

with the parameters as in eq. ( p.26|) . 

For |6| > 1, we must flip the fermion parity in the sector with negative 
phase shift a. However, as a result of the calculation, we found that the 
above equation ( p.33| ) is still valid in this case. Thus the partition function 
of the critical Ising model on the cylinder is given by eq. (|2.33|) for the entire 
range — oo < b < oo. Our result agrees with that by Henkel et. al. |T^], 
which is obtained in r-continuum limit. 



3 Identification of the boundary state by fold- 
ing 

In order to apply boundary conformal field theory to the present problem, 
we fold the system so that the defect line becomes the boundary of the 
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Defect Line 




Figure 2: The folding of the Ising model on a cylinder to a c = 1 theory on 
a strip. We fold at the defect line and also at the line on the opposite side. 
These lines correspond to the boundary in the folded system. 
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system [Q. The cylinder of circumference 1 is folded to a strip of width 
1/2, as shown in Fig. In addition to the boundary which corresponds to 
the defect hne, we have another boundary on the opposite side. The two 
boundaries are identical when there is no defect line, or equivalently when 
the periodic boundary condition is imposed on the Ising model (6=1). As a 
result of the folding, the degrees of freedom in the bulk are doubled. There 
are two different kinds of spin operators at each point, corresponding to the 
locations before the folding. We will call them two different layers of spins. 
Because of the doubling, we must consider a c = 1 conformal field theory, 
rather than c = 1/2. Known c = 1 conformal field theories are described by 
a free boson or its orbifold. 

Actually, a useful trick to study the Ising model on a plane is to regard 
two independent Ising models as a free boson theory [^, Let us denote 
the spins in each layer as cxi and a2- The product 0"i(j2 of the spin operators 
in each layer (at the same point) can be expressed as a simple bosonic op- 
erator cos 99 where 99 is the free boson field. Since the two Ising models are 
independent, correlation functions of this composite operator is always the 
square of the spin correlation in a single layer. Thus we can calculate spin 
correlation functions of the Ising model in terms of the free boson. In the 
present problem, the two Ising layers are coupled at the boundary though 
they are independent in the bulk. This makes the problem more difficult. 

The two-dimensional Ashkin- Teller model p5[ is defined by the classical 
Hamiltonian 

^ = Y1 {K[(Ji{,x)ai{y) + 0-2 (x) (72(1/)] + Lai{x)(7i{y)a2{x)a2{y)} , (3.1) 

where {x, y) runs over all nearest neighbor pairs on the square lattice. The 
doubled independent Ising model can be regarded as a decoupling point L = 
of the Ashkin- Teller model. Of course this is usually an unnecessary compli- 
cation when one focuses on the Ising model. However, in the present problem 
it seems necessary to consider the doubled Ising model as a special case of the 
Ashkin- Teller model. This will become manifest in Sec. where two-spin 
correlation functions are discussed. 

The critical Ashkin- Teller model is identified with a c = 1 conformal field 
theory. Precisely speaking, it is not a simple free boson but the Z2 orbifold 
of the free boson. (For a review, see Ref. |2^.) We take the normalization of 
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the free boson field so that the Lagrangian is given by 



C=Ud,^f. (3.2) 

ZTT 

The orbifold theory depends on a continuous parameter: the compactification 
radius r of the free boson (i.e. we make the identification ~ + 27rr). 
The decouphng point (doubled Ising model), which is relevant in our defect 
problem, corresponds to r = 1. The boundary condition (or the boundary 
state, see below) in the present problem can be identified in two (presumably 
equivalent) ways: in terms of Ashkin- Teller boundary states and in terms of 
the orbifold free boson. We first describe the former identification, which will 
be useful in the calculation of the spin correlation function. We will discuss 
the orbifold boundary state in the next section. 

A boundary of the system is described by a boundary state, if we exchange 
space and time so that the "(imaginary) time" is orthogonal to the boundary. 
In particular, when the boundary is conformally invariant, the boundary state 
\B) must satisfy the condition 

(Lf)-LL^j)|5) = 0, (3.3) 
where L*^-^) and L*^^) are the bulk Virasoro generators [^, |2^. A solution to 



this equation is given by the Ishibashi state 

|A)) = 5:|A,iV)®|A;iV) (3.4) 

N 

where | A, N) is a normalized A^-th descendant of the primary field of weight 
A, and the summation is over all descendants. In this paper, |A)) denotes 
the Ishibashi state constructed from the primary state with the weight A. 
A linear combination of Ishibashi states also satisfy conformal invariance. 
In general, the boundary state consists of several Ishibashi states. The cor- 
responding primary weights are included in bulk operator contents. While 
there are bulk operator with spin, only the spinless primaries (A = A) are 
relevant in the boundary state, as is seen from the structure of the Ishibashi 
states. 

The partition function and the operator content of the critical Ashkin- 
Teller model is studied by Yang ||30|. We summarize the spinless primaries 
at the decoupling point (doubled Ising model) in Table |1|. Hereafter we will 
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A = A 


Multiplicity 


Degeneracy 


Identification 


(n+l)2 
2 

(2n+l)2 

8 ^ 
(2n+l)2 

16 


1 
2 
1 

2 


[degenerate ] 
[non-degenerate ] 
[non-degenerate ] 
[non-degenerate ] 


h ■ h {n^ 0); ei • 62 (n = 1) 

ei • /2, -/"i -£2(71 = 0) 

(7i • (72 (n = 0) 

(7i • /2, -/"i • (^2 (ri = 0) 



Table 1: The spinless primary operators in the critical Ashkin-Teller model 
(at the decoupling point), n = 0,1,2,... and the "degeneracy" means the 
degeneracy in terms of the c — 1 Virasoro representation (not the multiplicity 
of the operators) . 



denote the critical Ashkin-Teller model at the decoupling point simply as 
"Ashkin-Teller model". 

A primary state |/) of the Ising model satisfies the condition 

4«-g|7) = (n > 0) (3.5) 

where Ljj^"^ jg ^ Virasoro generator of the Ising model. A primary state of 
the Ashkin-Teller model is defined by the condition 

^^Isingl ^ ^Ising2^) j^^.^ ^ q ^ ^g g^ 

Thus a tensor product of two Ising primaries is a primary (or a sum of 
primaries) of the Ashkin-Teller model, but the reverse is not always true. In 
fact, there is an infinite number of primaries in the Ashkin-Teller model, while 
only three primaries are present in the Ising model: the identity operator 
/, the energy density e and the spin operator a. For example, the second 
set of the above spinless primaries (A = A = (n + 1)^/2) is identified as 
follows. There are two primaries for each dimension. For n even (A = A = 
(2m -|- 1)^/2), the two Ashkin-Teller primaries are ei • I2 and Ii -62 {n = 0) 
or corresponding Ising descendants {n > 1). We will distinguish them as 
1^^^^^, 1)) and 2)). For n odd (A = A = 2m^), the two primaries 

are Ising descendants of Ji • I2 and ei ■62- We will distinguish corresponding 
Ishibashi states as |2m^,//)) = |2m^,l)) and |2m^,ee)) = |2m^,2)). 

Let us write the tensor product of the Ising Ishibashi states as Ashkin- 
Teller Ishibashi states. For example, we consider the Ising Ishibashi state 



14 



|/)) of the identity operator /. (Note that it is not a physically allowed 
boundary state of the Ising model without superposition with other Ishibashi 
states ^^.) The partition function on a strip of width 1/2 and length (3 with 



two boundaries specified by the Ishibashi state |/)) is given by the c = 1/2 
Virasoro character p9| 



xM = i^^WM + m)] (3.7) 

where q is defined as 

q = (3.8) 

Let us consider the Ashkin- Teller boundary state defined by the tensor 
product of the Ising Ishibashi state \l2))- The partition function of the 

Ashkin- Teller model on the strip with the corresponding boundary condition 
on both boundaries is given by the square of eq. ( p.7|) . Regarding the Ashkin- 
Teller model as a c = 1 conformal field theory, it should be possible to express 



it in terms of c = 1 Virasoro characters. The c = 1 Virasoro characters 
for the primary weight h are given by 

X,M) = ^, (3.9) 

ii h ^ n^/4. When h = n'^/A, the representation is degenerate and the 
character Xh is rather given by 

Xk{q) = ^- . (3.10) 

r]{q) 

The partition function on the strip is expressed as 

oo oo 

Xii^y = J2 X2n^{q) + J2 X{2n)2{q). (3.11) 

n=l n=0 

Thus the boundary state is identified in terms of Ashkin- Teller Ishibashi 

oo oo 

® 1/2)) = E |2^''^^)) + E 1(2^)'))- (3-12) 

n=l n=l 
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Similarly, we identify other tensor products as follows. 

oo oo 

|ei))®|62)) = Y.\'^n\ee)) + Y.\{2n+lf)) 

ra=l ra=0 
n=0 ^ 

\h))®\e.)) = t\^^^^^,2)) 

|ei))®|/2)) = El 9 .1)) 

n=0 



n=0,3{mod4) 

16 



I \\ ir \\ V- ,(2n+ 1)2 



n=0,3(mod4) 

16 



n=l,2{mod4) 

(2^ 

16 



n=l,2(mod4) 



The boundary states of the Ising model correspond to the free, up-spin 
and down-spin boundary conditions are given by PB| 

3.21) 
3.22) 

3.23) 



I/) = 

IT) = -^m + m + ^,w)). 
\i) = ^(i/)) + ie)))-^ia)). 



3.13) 
3.14) 
3.15) 
3.16) 
3.17) 

3.18) 

3.19) 

3.20) 



The 6 = defect in our model is described by the tensor product of the free 
boundary state of the Ising model as 

Iff) = (|/l))-kl)))®(|/2))-|62))) 

= j:t'^os^\'^,j)) + f:\n^)). (3.24) 
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Similarly, the fixed spin boundary states of the Ashkin- Teller model are given 
by 

ITT) = |Ti)®IT2) 

1 2 oo 2 1 oo 

^ j=l n=l ^ ^ n=0 

+ V2„^o' 8 ^^"^23/4^1.1 16 ,^)) (3-25) 

1 2 oo 2 1 oo 

lu) = ^EElT,i-)) + ^EI^')) 

^ j=l n=l ^ ^ n=0 

ITi) - ^EEIy,j-)) + ^EI-^))-4^EI^^^)) 

j=ln=l n=0 V n=0 

^ E(-iyEl^^^^,j)) (3.27) 



UT) = ^EEIy,j)) + ^EI-^))-4ei^^^)) 

^ j=l n=l ^ ^ n=0 V ^ „=o ° 

' E(-l)^'El^^^^.i)) (3-28) 



Furthermore, there are tensor products of free and fixed boundary states: 
IT/> = \f.\'^,A)) + \(-lT\n')) 

l^.n^ ... 1 



IT/) = 2Y.\^^^)) + 2^-mn')) 



n=l 



(3.30) 



n=0 
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where 



l/T) 


1 oo 

= ^E(-i) 

^ n=l 


' 2 




-| oo 

1 V 

' 91/4 

^ n=0 




l/i) 


-1 00 

= ^E(-i) 

^ n=l 


' 2 




00 

9I/4 E 
^ n=0 


(-1) 




1 

2' ^ V2 









n(n 



±llA2n+lf 



16 



.2)) 



(3.31) 



^)) + ^(-ir 

^(n+i) (2n + 1)" 



, n 



16 



,2) 



2)) 



(3.32) 



(3.33) 



It is natural to think that the combination 

I TT) 



2 °° .n' 



ii) = E E I'i^j)) + v^E 1^^^^)) + E I-')) (3-34) 

j=ln=l ^ n=0 ° n=0 



corresponds to the infinitely ferromagnetic defect b ^ 00 together with the 
anisotropic limit K2 0, namely a = —j. (The overall normalization is 
determined from the partition function. See below.) Similarly, 



2 - ,n2 



Ti) + 1 it) = E E I'i^j)) - v^E 1^^^^)) + E 1^')) (3-35) 

j=ln=l ^ n=0 ° n=0 

corresponds to the infinitely antiferromagnetic defect b —00 with K2 0. 

Thus we have obtained the boundary state for three special types of defect 
line. Now we make the following assumption about the boundary state for 
general strength of the defect, generalizing ( p.24| ),( |373^ and ( |3.35| ). 



Assumption 

The boundary state \B) is given by 

00 00 2 

\B{vo)) = E \n')) + V2j:cos{n^o)\^,S)). 



n=Q 



Here ipo is given by 



^0 



n=l 

TT a 
4^2' 



(3.36) 
(3.37) 
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using a defined in ( ^.Idj ). \n'^/8,S)) is defined as 




We note that the boundary state contains only the combination |n^/2, S)) 
of Ishibashi states, which is symmetric with respect to the interchange of 
the two Ising layers. This is presumably a consequence of the reflection 
symmetry about the defect line in the unfolded picture. For the boundary 
corresponding to the defect line, we take a = a{b). For the other boundary, 
we take ipo = vr/4 because there is no defect {b = 1). We summarize the 
correspondence between the defect strength and the parameter (po in Table ^ 





b 


Condition 


Description 





OO 


K2^0 


Infinitely ferromagnetic limit 


7r/4 


1 


(none) 


No defect / periodic 


tt/2 





(none) 


Free boundary 


37r/4 


-1 


(none) 


Antiferromagnetic defect / antiperiodic 


TT 


— OO 


K2^0 


Infinitely antiferromagnetic limit 



Table 2: The value of ipo for several strengths of the defect 



We can calculate the partition function on the strip from the boundary 
state. Using the transfer matrix H^^^ in the direction orthogonal to the 
boundaries. 



Z 




OO CO 



(g) + 2 ^ COS {nipo) COS [nir/ 4)x„2/8 (q) 



n=0 n=l 



Z{^o--r) + Z{ipo + -) 



(3.39) 



where 



z{0 
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m) 

= ^^"''^'^"[^3(e-^«^ q) + Ue-^^^ ?)] (3-40) 

This actually agrees with the partition function of the model obtained in ( |2.33| ) 
We note that, while Henkel et al. ||12|, |13[ identified the partition function 
with a c = 1 Virasoro algebra empirically, it is a natural consequence from 
our "folding" approach. They further studied the system with many parallel 
defect lines [|^. For example, for the system with two defects, they found 
that the spectum can be described in terms of a Virasoro algebra if the 
location of the defect lines is commensurate, namely the distance between 
two defects is an integral multiple of 1/n of the system size, where n is an 
integer. The central charge in this case is given by c = n when n is odd 
and c = n/2 if n is even. This can be also naturally understood in terms of 
multiple folding, i.e. we fold the system many times so that we have 2c layers 
of the Ising model, and defect lines are placed only at the boundaries. In 
principle the boundary conformal field theory for the corresponding central 
charge could be useful for analysis of the commensurate multiple defects. 
Such an analysis is however beyond the scope of the present paper. 



4 Z2 orbifold of free boson and the boundary 
state 

We turn to the second identification of the boundary state in terms of the free 
boson. This is useful in calculation of some correlation functions, as we will 
see in the Section ^. Moreover, it gives a more transparent understanding of 
boundary states and also enables us to find a new universality class of defect 
lines. 

4.1 Boundary states of the free boson 

In order to make our paper self-contained, here we give a brief review on the 
boundary states of the free boson. We consider the free boson (before the 
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orbifolding) defined by the Lagrangian density 

^ = ^{d.^r (4.1) 

where ip = <f{cr,t) and < a < (3. (f is compactified with the compactifica- 
tion radius r, namely ~ y9 + 27rr. From the Lagrangian, we can derive the 
equation of motion 

9V = (4.2) 
and the canonical commutation relation 

b(a,t),nK,t)] =25(a-a') (4.3) 

where 11 = 0/vr. We impose the periodic boundary condition in "space" 
(parallel to the boundary direction) ip{<J,t) ~ ip{a + P,t). 

We can determine the mode expansion of the boson field from eqs. ( |4.2| ) 
and ( [4.3|) as follows: 

ip{a,t) = 

Stt tt „ 

X + —rwa + —pt 

P P 



where w is an integer (winding number) allowed by the angular nature of the 
boson field. The operators satisfies the commutation relations 

= i (4.5) 

[ttn, a}ra\ = [«n, = ^nm (4.6) 

and the other commutators vanish. Since the constant mode x is also com- 
pactified as X ~ X + 27ir, the conjugate momentum p is quantized to an 
integer multiple of 1/r. 

The boson field if can be decomposed into chiral components cpi and (pn 

as 

^l{x^) = - + -{rw + -)x+ 
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1 °° a 



(4.7) 



(4.8) 



where = t ±a. 

We make a mode expansion of the energy-momentum tensor T{x^ 
^{d+cpLT and f(x-) = ^{d^cpRf as 



T{x^) 



27r 

2% 



m=— oo 
oo 



and L^, which are (bulk) Virasoro generator, are given by 



n=—oo 
oo 



1 



n=— oo 



where a™ and a™ are defined as 



ar, 



—i\fnan {n > 0) 

(rw + |) {n = 0) 
. ^A/nal^ (n < 0) 



and :: denotes the normal ordering 



iy/ncLn {n > 0) 
i-rw + D (n = 0) 
(n < 0) 



amO-n {n > m) 



(4.9) 
(4.10) 



(4.11) 
(4.12) 



(4.13) 



(4.14) 



(4.15) 
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The Hamiltonian for the free boson with the periodic boundary condition in 
the spatial direction of length /3 is given by 



rr(P) _ — 



27r 



(4.16) 



A conformally invariant boundary state must satisfy the Ishibashi condi- 
tion (|3.3|). In the present case, this condition is satisfied, if the boundary 
state \X) satisfies 

(a„±a_„)|X) = (4.17) 

for any integer m. (The sign ± should be common for all m.) Let us first 
consider 

(«„ + a_^)|X) = 0. (4.18) 
The conditions for m 7^ are satisfied by 

00 

|X)=exp[-5:atat]|vac), (4.19) 

n=l 

where |vac) is an oscillator vacuum which satisfies a„|vac) = a„|vac) = 0. An 
oscillator vacuum is characterized by two integers: w and k = rp. (Recall 
that p is quantized due to the compactification.) Hereafter \{w,k)) denotes 
the oscillator vacuum with zero-mode parameters w and k. The remaining 
condition for m = requires w = 0. Thus a boundary state which satisfies 
eq. ( |4.18| ) is given by 



Y^C,exp[-J2al~ai]mk)), (4.20) 



n=l 



where Ck are constants. Besides the Ishibashi condition, there is a consis- 
tency condition found by Cardy We can calculate the partition function 



on a strip with two boundary conditions. If we exchange the space and time 
and regard the direction parallel to the boundary as time coordinate, the 
partition function should be expressed as a sum of Virasoro characters with 
integer coefficients. Namely, the partition function with the boundary con- 
ditions A and B on the two sides should take the form 

Z = j:<BXhiq) (4.21) 

h 
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where are non-negative integers and Xh{(l) is the Virasoro character for 
the primary weight h. If the two boundary conditions are the same, the par- 
tition function shows the boundary operator content for the boundary. In 
particular, there should be exactly one dimension- zero character correspond- 
ing to the identity operator: = 1. These requirement comes from the 
radial quantization on a half-plane and the conformal mapping to the strip. 
A solution to the consistency condition is given by 

Ck = ^e-*'^^"/", (4.22) 
V2r 

where cpo is a constant. It gives a one-parameter family of boundary states 

1 oo oo 

\D{vo)) = ^ E e-^'^^/^exp[-j:al~aimk)), (4.23) 

Below we show that this state satisfies the consistency condition. Let us 
assume the width of the strip is 1/2 (and the length is P). The partition 
function of the strip for the above boundary condition with parameters ipo 
and (fQ at two boundaries is 

Z,(A^o) = {Di^o)\e'<'/'\Dicp',)) (4.24) 



1 1 



2r 77(g) 



where q is defined in eq. ( p.8| ) and Aipo = ipo — ipQ. Here l/r]{q) is given by the 
summation over the oscillator states and 7^3 part comes from the summation 
over zero-mode quantum number k. We express the above partition function 



as a function of q defined in eq. ( p.26| ). This is achieved by a modular 



transformation of rj and function. The result is 

Z,(Av.o) = ^g-(A^oA)^^3(e-4.A^o/3^ ^4.^)^ (4.25) 
7]{q) 

This actually is a sum of c = 1 Virasoro characters with non-negative integer 
coefficients. In particular, it gives the boundary operator content when the 
two boundary states are the same {Aipo = 0). In this case, the partition 
function contains one dimension-zero character. Thus Cardy's consistency 
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conditions are satisfied. The boundary state ( ^4.23| ) has a simple physical 
meaning: we see that 

<^(cr,t = 0)\D{ipo)) = Vo\D{ipo)), (4.26) 

namely the boson field takes a fixed value ifQ at the boundary. Thus it 
corresponds to the Dirichlet boundary condition ip = ipo. 
We can discuss the other possibility 

(a„-«_„)|X) = (4.27) 

in a similar manner. In this case, conditions for m ^ implies 

oo 

|X)=exp[+^44]|vac), (4.28) 

n=l 

and that for m = requires p = 0. A solution to Cardy's consistency 
condition is given by 

oo oo 

|A^(^o)) = v^ E e-2-^°exp[+E44]l(«^,0)). (4.29) 

ui=— oo n=l 

Again it has a simple physical interpretation: von Neumann boundary condi- 
tion (0 = const.), or equivalently, Dirichlet boundary condition on the dual 
field (p = ifL — ifji = (pQ. Since the von Neumann boundary condition is the 
Dirichlet boundary condition on the dual field, the amphtude between two 
von Neumann boundary states are given by replacing the compactification 
radius by its dual l/(2r), and the parameters ipo and ip'^ by (pQ and ip^ in 
eq. ( [4.25| ). Thus the von Neumann boundary state is mutually consistent 
with a von Neumann boundary state with any value of the parameter. We 
note that, although the Dirichlet and Neumann boundary states contain a 
continuous parameter, the operator content of the theory does not depend 
on the parameter, as is seen from the amplitude with the same boundary 
state. 

The amplitude between the Dirichlet and the von Neumann boundary 
states is given as 

(p\ 1 oo -| 

(D(^o)ie--^^/^iiv(^o)) = Tfr^/^nrr^ 
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(4.30) 



This is again a sum of c = 1 Virasoro characters with non-negative integer 
coefficients; they are mutually compatible. We note that this Dirichlet-von 
Neumann amplitude depends neither on ip^ nor (p^, because only the zero 
winding number sector built on |(0,0)) contributes to the amplitude. 



4.2 Boundary states of the orbifold 

Now we discuss the boundary states of the Z2 orbifold, which has a direct 
relevance to our problem. Related discussion in string theory can be found in 
Refs. [p2| , |33| , |3^ . The Z2 orbifold of the free boson is constructed from the 



free boson, identifying ip ~ —tp. Thus a boundary state must be invariant 
under the transformation 

G:^^-^. (4.31) 

The simplest way to satisfy this requirement is to symmetrize the free bo- 
son boundary state. Such a boundary state constructed from the Dirichlet 
boundary state of the free boson is 

\Do{^,)) = ^m^,)) + \D{-^,))1 (4.32) 

where < < vrr. The overall factor is determined from the consistency 
condition. The partition function of the orbifold theory for the Dirichlet 
condition with parameters (p^ and ip'^ at the two boundaries is given by 

{Do^ipoW'^''' '^\Do{v'^)) = Zr{ip^ - ip'^) + Zr{ipo + v'o), (4.33) 



where Z^ is defined in ( |4.25| ). Thus the orbifold Dirichlet boundary state (|4.32| ) 



is mutually consistent boundary state for any < yjo < ^rr. The operator 
content for the boundary state ( [4.32| ) is obtained from the above partition 
function with setting ip'^ = ip^. We note that, it contains Z^{2ipQ) which 
depends on the continuous parameter ifQ. Thus, the boundary operator con- 
tent of the orbiford theory depends on the boundary value paramter yjQ- On 
the other hand, the boundary operator content for the Dirichlet boundary 
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state ( |4.23| ) of the (unorbifolded) free boson does not depend on the bound- 
ary value ifo- This is a consequence of the U{1) symmetry (f ^ (f + const, of 
the free boson. The orbifolding breaks this symmetry and thus the spectrum 
of the boundary operators can depend on ipQ. 

Similarly, we can consider the orbifold version of the von Neumann bound- 
ary state: 

\Noi^o) = ^ [|iV(^o) + liV(-^o))] , (4.34) 



where |iV((^o)) is given in ([4.29|) and < (fo < n / (2r). It can be also shown 
that it is mutually consistent boundary state, and also consistent with the 
orbifold Dirichlet boundary state (|4.32|) . From the analysis of the partition 
function, the operator content of the state ( [4.34| ) also depends on (pQ. 

The endpoints (po = 0,7rr and (po = 0,7r/(2r) are the fixed point of the 
transformation ( (4.311 ). At those endpoints, the orbifold boundary state ( [4.32| ) 
or (|4.34| ) contains two dimension- zero boundary operators and therefore does 
not satisfy Cardy's consisntency condition. On the other hand, another kind 
of Dirichlet /von Neumann boundary state is possible at these endpoints. 

Namely, due to the orbifold condition ip ~ —ip, the antiperiodic boundary 
condition in terms of free boson 



(p{a + P,t) ^ia,t) 



(4.35) 



is allowed under the "periodic" boundary condition of the orbifold. Let us 
call the subspace with the above boundary condition the twisted sector. The 
mode expansion of the free boson with the antiperiodic boundary condition 
is given by 



^{a, t) 



^ n=\ 



^ n=l 



^ A Un-l){a+t)^-f 



n — 



i(n-i)(cr-t) 



n 



p + 



hi 



n — 



(4.36) 



where bnfilifin and b}^ are the boson creation/annihilation operators obeying 
the commutation relations similar to (14.61). The Hamiltonian for the free 
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boson with the antiperiodic boundary condition in the spatial direction of 
length 13 is given by 



Tr{A) _ — 



27r 



(4.37) 



The essential point in the twisted sector is that x can only take the fixed 
point values or vrr due to the the antiperiodic boundary condition ( |4.35|) . 
Thus there are two independent oscillator vacua |0)t and iTrr)^- corresponding 
to a; = 0, 7rr, in the twisted sector. From them we can construct the Dirichlet 
boundary state in the twisted sector. Those twisted Dirichlet boundary states 
are given by 

i^v^u;./ ^ "-i''"'"|¥'o)t, (4.38) 

where (yJo = or vrr only. 

The von Neumann boundary states in the twisted sector are less clear. 
Since if is not fixed in the von Neumann boundary condition, we assume that 
they are constructed on the oscillator vacua |0)T±|7rr)T. We will later confirm 
that this ansatz actually leads to mutually consistent boundary states. Then 
the twisted von Neumann boundary states are given by 

|Ar(0)r) = e+Sn=i''"''"-^(|0)T+kr)T), (4.39) 

\N{^)t) = e-^T.ZAiil=(\Q)^-\^r)T). (4.40) 

While these four twisted boundary states satisfy the Ishibashi condition, 
they are not consistent boundary states by themselves. The amplitudes 
among them are 

{D{^.)T\e-<'"\D{^.)T) = g"/^^n^4n7^ 

„=i 1 - 



2r/(l) 

(iV(^o)T|e-<'/'|iV(^o)T) 



(4.41) 
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= ± 



- ±^^^q"''m-"\q'l (4.42) 
v2?7(g) 

where = or 7rr, (^0 = or tt/ (2r) and ± depends on the value of (/?o and 
(pQ. The other amphtudes vanish: 

(D(0)T|e-<'/'|i^(7rr)T) = (A^(0)T|e-<'/'|7V(|:)T) 

= 0. (4.43) 

In order to construct consistent boundary states, we must combine the 
twisted and untwisted sectors with appropriate coefficients. These endpoint 
Dirichlet/Neumann boundary states are 

\Do{v^)±) = 2-^/2|^(^o))±2-^/1/^(vpo)t), (4.44) 
\No{^o)±) = 2-'/'\N{ipo))±2-'/'\N{ipo)T), (4.45) 

where (po and (^0 again take only the endpoint values. The amplitude among 
those endpoint Dirichlet boundary states are 



2 " ' 277(g) 

= ^ E K"' + (-1)" 



Ml) „±-, 

00 

Yl X2r2n2 + X(2n)2 , (4.46) 



n=— oo 

oo 



n=l n=0 



(L>o(¥'o)±|e-<'/'|£'o(¥'o)T) = Iz^iO) - -^Mq') 



2 ' 277(g) 
-1 00 

00 00 

1] X2r2n2 + X] X{2n+1)2 , (4.47) 



n=l n=0 



29 



_tt(A) 10 1 

{Do{vo)vi\e ''\Do[T^r ~^q)V2) = -^r(vrr) 

oo r2(2n-l)2 

= E / N . (4-48) 

n=l ^(9) ^ ^ 

where z/i 2 take either + or — independently. The amphtude among the 
endpoint Neumann are given by replacing the radius r by its dual l/(2r). The 
amplitudes between the endpoint Dirichlet and Neumann boundary states are 
given by 

Ar]{q) 27]{q) 

00 

= EXi(2n-l±i)2, (4.49) 
n=0 

where ± in the last line depends on the signs z/i 2 and values of ipo and 
ipQ. Thus these eight discrete boundary states mutually consistent. More- 
over, they are consistent with the continuous Dirichlet ( [4.32| ) and von Neu- 
mann ( [4.34| ) boundary states on the orbifold. This can be easily seen because 
only the untwisted sector in ( ^.44| ) and ( ^.45] ) contribute to the mutual am- 



plitude with the continuous ones. For example, 

(Do(0) ± le"^/' ^'\Do{^o)) = Zriipo). (4.50) 

In this way, we have found two continous families of boundary states 
and eight discrete ones and they are all mutually consistent boundary states. 
While our construction does not exclude the presence of other boundary 
states, we conjecture that the above exhaust all the possible orbifold bound- 
ary states for generic values of the compactification radius r. An analogous 



statement has been proven for a non-orbifold boson [36|. (For special val- 
ues of r, another family of boundary states is constructed for non-orbifold 
boson Similar construction may apply to the orbifold boson.) 



Finally, we note that although the above orbifold boundary states are mu- 
tually complatible, they are incompatible with non-orbifold boson boundary 
states. For example, the coefficient of the boson Dirichlet boundary state is 
fixed as ( |4.23| ), by requiring the compatibility with itself. Thus the partition 
function with the boundary conditions (|4.32| ) and ( [4.23| ) at the opposite sides 
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is given by 

--^[Zr{(fo-f'o)+Zr{fo + f'o)]. (4.51) 

This contains non-integer (irrational) coefficient when expressed as a sum of 
c = 1 Virasoro characters; the two boundary conditions are incompatible. 
This incompatibility is presumably related to the fact that the sets of local 
operators are different for the free boson theory and its orbifold (Ashkin- 
Teller model). 



4.3 Applications to the defect line problem 

Now we turn to our original problem of the Ising model with a defect line. 
As explained, the bulk theory in the present problem is the special r = 1 
point of the orbifold, after the folding. Thus a universality class of the defect 
line should correspond to an orbifold boundary state with r = 1. Hereafter 
we fix the compactification radius as r = 1. 

From the general discussion of the orbifold boundary states, we have two 
continuous family of boundary states and eight discrete ones, all of which 
are mutually consistent. The continuous Dirichlet boundary state ( [4.32| ) is 
identified with the boundary state (|3.36|) for the defect line in ( p.lD , using 
the same (po defined in ( |3.37| ). In fact, the partition functions constructed 
from ( p.36|) and (|4.32| ) agrees completely, observing that Z{Alpq) in eq. ( p.40| ) 



is identical to Zr{A{po) in eq. (^4.25|) for r = 1 



Comparing ( |3.36| ) and (f4.32|) , we can identify the symmetric combination 



of bosonic oscillator vacua as 



[1(0, A:)) + 1(0, -A:))] = 1^,^), (4.52) 



where the right hand side is defined as eq. (|3.38|) in the Ashkin- Teller classi- 



fication. The other symmetric combination is also identified with eq. ( p.33| ) 

as 

1 

-^[|(n,0)) + |(-n,0))] = |y,A). (4.53) 

Other 3x3 = 9 universality classes can be readily obtained by cutting 
the system at the defect line, and imposing one of the three possible Ising 
boundary universality class on the two sides. As discussed in the previous 
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section, the corresponding boundary states are given by the tensor product 
of the Ising boundary states. One of them, |//) is actually the special case 
= vr/2 of the Dirichlet boundary state ( |3.36D or ( |4.32| ). The remaining 8 
universality classes nicely fit with the 8 discrete endpoint boundary states of 
the orbifolds. Actually we can identify them as following: 



ITT) = 


l^o(0)+) 


(4.54) 


ITi) = 




(4.55) 


UT) = 


\Do{^)+) 


(4.56) 


lii) = 


l^o(O)-) 


(4.57) 


T/) = 


|iVo(0)+) 


(4.58) 


if) = 


WoiO)-) 


(4.59) 


l/T) = 


|iVo(|)+) 


(4.60) 


l/i) = 




(4.61) 



We note that the sign flip of the coefficient of the twisted sector corresponds 
to an overall Ising spin flip. 

These are equivalent to the following identification between the twisted 
orbifold boundary state and the Ashkin- Teller Ishibashi states: 

(2n+l)2 (2n + l)2 ^ 



\D{0)t) 


V2 „=o \ 






\D{n)T) 


1 ^ / 

^2 n=0 \ 


|A^(0)t) 


oo 

= E(-ir 




n=0 




oo 

= E(-ir 

n=0 



16 

(2n + 1)2 



16 



1)) + 
1))- 



16 

(2n + 1)2 



16 



,2)) 



(2n + l)2 
16 

(2n + l)2 
16 



1) ) 

2) ) 



(4.62) 
(4.63) 
(4.64) 
(4.65) 



These correspondences are verified from the calculation of the partition func- 
tions: 

oo (2n + l)2 - , 1^2 



v2 „=o 



E((- 

n=0 

(2n+ 1) 
16 



16 



■|(-1) 



L) |^-(Lo+Zo)/2| (2n + l)\ ^ 



1 



16 

■'("+!) . {Lq+Lq) (2n + 1)^ 

2 q 2 \i i- 

^ ' 16 
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which should be compared with ( |4.41| ) and ( [4.42| ). 



Thus all the known defect universality classes are understood in terms of 
orbifold boundary states. On the other hand, the defect line corresponding 
to the continuous von Neumann boundary state ( |4.34]) has not been yet 



identified. This actually gives a previously unknown universality class of 
defect lines. The nature of this defect universality class will be discussed in 
Section 

5 Correlation functions for the continuous Dirich- 
let universality class 

Here we discuss several correlation functions in the presence of the defect 
line that is identified with the continuous Dirichlet boundary state ([4.32|) 
or (|3.36|) in our approach. 

5.1 Spectrum of boundary operators 



As discussed in Section [4.2| , the complete spectrum of boundary operators 
can be determined from the partition function of the strip with the same 



boundary condition at the both sides ||2^. In the case of the continuous 
Dirichlet universality class ( |4.32| ), it is obtained by putting (p'^ = (po and 
r = 1 in eq. ([4.33|) as 

Z = Z(0) + Z(2(^o), (5.1) 
where Z(A(^o) is defined in ( p.40| ). By expressing this result in terms of c = 1 



Virasoro characters (3^), we obtain the complete set of boundary primary 
operators. 

The boundary operators consist of two sets: the one corresponds to Z{0), 
which is independent of the boundary condition (defect strength), and the 
other corresponds to Z{2{po) which depends on the boundary condition. As 
we discuss later, the former is identified with the boundary operators asso- 
ciated with the bosonic operators, and the latter is the boundary operators 
associated with twist (spin) operators. The surface critical exponent of the 
operator is defined by the correlation function of the two operators located far 
apart along the boundary (the defect line in our case). This surface critical 
exponent is given by the dimension of the corresponding boundary operator. 
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Here we give the complete spectrum of the boundary primary operators in 
the "moving sector" Z{2ipo): 

A. = 2n^ + l!^ + 2^, (5.2) 

where n is an arbitrary integer. (The above is the complete set for generic 
!fQ, where no degenerate representation of c = 1 Virasoro algebra occurs.) 

Let us assume that the surface exponent of the Ising spin is given by 
the boundary operator with the smallest dimension in the above "moving" 
sector. Then we find that the exponent depends continuously on the defect 
strength as 

2A, = minf4(^)^4(^^)^). (5.3) 

This continuously changing surface exponent reproduces the previous result 
by Bariev and McCoy and Perk |^. (Note that their lattice model is dual 
to our model, and surface exponent changes under the duality transformation 
as discussed in Ref. p. In Section we argue they belong to the same 
universality class with different parameters.) 

Our formalism gives a complete set of boundary operators in addition to 
their result. Note that the surface exponent is unchanged under (po ^ n — ipo 
which corresponds to the sign reversal of the defect b — *■ —b. 

We note that, for = or vr there is a dimension zero operator other than 
the identity. Thus Cardy's consistency condition does not hold in these cases. 
The presence of the dimension zero operator is related to the infinitely long 
range correlation of the spin operators along the boundary. This somewhat 
unphysical situation is only achieved as a limit of anisotropic and strong 
defect limit. (However, it corresponds to 6 = in dual quantum model ( |5.46| ), 
which is a well-defined quantum model.) 



5.2 Two-spin correlation functions 

The conformal field techniques enable us to calculate also the exact two-point 
correlation functions at general location for arbitrary strength of the defect. 
Let us consider the two-spin correlation function. By folding, the correlation 
of two spins located on the same or opposite side of the boundary is mapped 
to the correlation function (cticti) or (o"i(T2), respectively. Let us take the 
boundary as the real axis of the complex plane. In boundary conformal field 
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theory, a basic technique is to analytically continue the energy- momentum 
tensor beyond the boundary. Then the two-point correlation function of a 
non-chiral operator is mapped to a four-point correlation function of chiral 
operators as 



{aj{zi,zi)ak{z2,Z2)) = {aj{zi)aUzl)akiz2)a',,{z2)) , 



(5.4) 



where z^ and Z2 are complex conjugates of zi and Z2, and (yj{zi) is a chiral 
operator which should be distinguished from the non-chiral aj{zi,zi). This 
may be regarded as a generalization of the mirror image method. 

In the present problem, the non-chiral spin operator aj {z, z) has the di- 
mension (1/16, 1/16) and all the four operators in (|5.4| ) have the dimension 
1/16. Thus the correlation function should be constructed from the c = 1 
conformal blocks of four spin operators of dimension 1/16. These conformal 
blocks are obtained by Zamolodchikov [^] as 



/(A,-,l,x) = (16v^)^^3i(n) 



(5.5) 



where A is the dimension of the intermediate primary state, x is the cross- 
ratio 

' Zi- Z2 



X 



and u is determined from x by 



X 



Zi — Z2 



(5.6) 



(5.7) 



(Note that our definition of the theta functions is different from that used in 
Ref. 1^.) Any real dimension A of the intermediate state is allowed in this 
c = 1 theory. 

The correlation function can be determined if the intermediate states 
and the corresponding OPE coefficients are identified. These OPE coeffi- 
cients depend on the boundary condition (defect strength) and in general 
different from those in the bulk four-point function. Nevertheless, Cardy 
and Lewellen |jl9[ presented a method to determine them from the bulk four- 
point function and the boundary state. We use this idea to determine the 
boundary correlation function (| 
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The essential point of their argument is that the OPE of the two bulk 
operators is independent of the boundary condition. If we use the bulk OPE, 
the boundary correlation function is in principle obtained from the one-point 
function in the presence of the boundary: 

{aj{zi, zi)akiz2, Z2)) ~ C<Tja,(A,A)(0A.A + (descendants)), (5.8) 

A, A 

where 0a, A is a primary field with dimension (A, A), Co.^.o-j.(a,a) are the bulk 
OPE coefficients and the coordinate dependent factor is omitted. The one- 
point boundary correlation function of the primary field 0a, A is given in 
terms of the boundary state as 

(0a,a(i/)) = 




(5.9) 

where y is the distance from the boundary. The boundary state \B) consists 
of Ishibashi states. In general, only the primary states \4>aa) with A = A 
can have non- vanishing A^. Once the contribution from an intermediate 
primary state is identified in this way, contributions from its descendants are 
determined from the conformal invariance. The total contribution from a 
conformal family should be the conformal block function multiplied by the 
constant A'^, as required from eq. ( ^.4| ). 

Zamolodchikov also constructed an explicit expression of the bulk 
correlation function of four spin operators of the Ashkin- Teller model in the 
same layer. In the doubled Ising case, it reads 

(0-1(2:1, ^1)0-1(2:2, ^2)0-1(^3, ^3)0^1(2:4, Z4)) = 

which is equivalent to the well-known four-point function of the Ising model. 
From this correlation function, we can read off the dimensions of the inter- 
mediate primaries as 

X^ . ,9 ,9^ A^^ + (2/+l)^ 

(A,A) = (2A;^2/2) C ^ , o )• (5-11) 
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In the calculation of the boundary correlation functions, we need only the 
spinless primaries 

A = A = ^ (5.12) 

and the corresponding OPE coefficient. 

To calculate the correlation function of the spins in the opposite side of 
the defect line {{aia2) after the folding), we need the OPE of (Xi and (T2. 
This can be read off from the bulk four-point function (o"icr2(Ticr2). Since the 
two layers are not coupled in the bulk, this is simply a product of two-point 
functions (cticti) {cr20'2) ■ Still it can also be written in terms of the c = 1 
conformal block as 

((Ji(zi, Zi)a2{z2, ^2)0"! (23, ^3)0^2(^4, Z4)) = 

Z14Z32 2 ^ W2fc-l)V8 W2/-1)V8 

TTT 2^ V M VM • 5.13 

Thus the spinless primaries appear in the OPE of cricr2 are 

A^A^i?^. (5.14) 

There is some ambiguity in determination of the OPE coefficient from 
the bulk correlation functions. The sign of the OPE coefficient cannot be de- 
termined from the bulk correlation where the square of a coefficient appears. 
The sign actually depends on the definition of the intermediate field, but the 
boundary correlation function depends only on the product A^C^cr^, which 
is independent of the sign convention of the field 0. In the following, we fix 
the sign convention so that > 0. Anyway we cannot determine the sign 
of A^Co-cr^ from those bulk correlation functions. 

Thus we assume a specific form of the OPE coefficients that is consistent 
with the bulk correlation functions, and construct the boundary two-point 
function. Then we will check the validity of our assumption by comparing 
the result with the known cases. 

We assume all the OPE coefficients are positive. The spinless operators 
appearing in the OPE of ai and (J2 have dimensions as shown in eq. ( p.l4| ) 
and are identified with the third set of the Ashkin- Teller primaries listed in 
Table ^ The bulk correlation function (|5.13| ) uniquely determines the OPE 
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coefficients as 

/ ]^ ^ (2n-l)2/8 

Cai<T2(2n-l)2/8 = v2 ( — ) (5.15) 



16/ 



where (2n — 1)^/8 is the primary weight of the intermediate operator. The 
spinless operators in the OPE of cTj and CTj (j = 1 or 2) have dimen- 
sions ( p.l2| ). They are naturally identified with = in the first set and 
the whole second set. This is consistent with the Ising OPE: aa ~ [J] + [e]. 
The second set with the primary weight has multiplicity two for each 
dimension, as we discussed in the last section. We determine which opera- 
tor appears in the OPE, from the Ising OPE. For odd n = 2k + 1 (weight 
(2A;+ 1)^/2), they are identified as the descendants of ei®/2 and Ii®e2- Thus 
the OPE of (Ji and cxi produces the former, and that of (T2 and (T2 produces 
the latter. For even n = 2k (weig ht 2fc2), they are identified as the descen- 
dants of Ji I2 and ei ® 62- From the Ising OPE, we see that both OPE of 
(Tj and aj (j = 1, 2) produces the former. Since the boundary state ( p.36| ) 
consists only the symmetric combination of corresponding Ishibashi states as 
in ( |3.38| ), we consider the OPE coefficient for the symmetrized operator. It 
is determined from the bulk correlation function ( 5.101) as 



Ca,a,(ny2,S) = V2i^ — j (5.16) 

for j = 1, 2 and n > 1. The coefficient for identity operator is 

a,.,o = 1 (5.17) 

for j = 1, 2. 

The coefficients is read off from ( p.36| ) as 



AnVs,s)_\ 1 {n = 0) 

v^cos(n^o) {n>l) ^^'^^^ 



Therefore we obtain the boundary correlation functions as 

/ \ 1/8 -j^ 00 2 1 
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1/8 -I 

t93(e'^^",v^), (5.19) 



\Ay1y2x J ^s{u 

(0-1(2:1, ^1)0-2(2:2, ^2)) 

1 ,(2n+l)V8^ ..(2n+l)^ 1 

: :^7::tz^^b c,,,,(2„+i)2/8/( — - — ,T^,i>a;) 



^92(e'''^°,v^), (5.20) 



1/8 



where yi and ?/2 are the imaginary parts of zi and 2:2, a; is the cross-ratio 
and M is a function of x determined by eq. (|5.7|). This generahzes Brown's 
result [0, which is the first order perturbation in the defect strength. 

Let us check the result, in order to justify our assumptions about the 
OPE coefficients. 

No defect (b=l) 

(^1^1) = 1 ~lTun — ^ M/8 Q ] ^ ^3(e'^/^ Vu) 

\zi - 22^1^ (1 - xyl^ ifslu) 



1 



(0-1^^2) 



\Zl-Z2\'/' 

1/8 

^^?2(e^^/^v^) 



(5.21) 



1^1 -^2*P/^ 



x{l — x) 



1 1 



1 



^2(e*"/^v^) 



|Zl-^2*|V4 



(5.22) 



Thus we recover the known result (bulk two-spin correlation) for both (aicri) 
and (cri(T2). 



39 



Free boundary {b = 0) 



(o-io-i) 



1 



1 1 



(0"l0"2) 



1 
1 



=^4(v^) 



-1/4 _ r^l/A 



X' 



\Z\ — Z2 



11/4 



xil — X) 







(5.23) 



(5.24) 



The former agrees with the boundary two-point correlation of the Ising model 
with the free boundary condition p|. The latter result is expected since the 
two layers are completely decoupled at this point. 

Antiferromagnetic defect (6 = —1) 



(cTlCTi) 



(0-lC^2) 



1 



1 



Zl-Z2|l/4(l-x)V8 §^{u) 



^^3(e^=^-/^v^) 



ki-^2|^/^ 
1 

1 

1 



x{l — x) 

1 



1/8 



^2(e*'"/',v^) 



1^1 ~ ^2 I 



^92(m)^4(m) 



(5.25) 



(5.26) 



The former is the same as the no defect case (6 = 0), while the latter changes 
sign. This is consistent with the original lattice model, where one can reverse 
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spin operators 




boundary 

Figure 3: The location of spin operators used in the graph of correlation 
functions (Fig. ^). We fix the distance of spin operators from the boundary 
to 1, and leave the horizontal distance r between them as a variable. 
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(1)) CV)rr(4a.ti()ii function of two spins in diffbront la\"ors 
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(po = O.Itt ; 
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Figure 4: The two-spin correlation for various strength of the defect. We 
show the result on (a) (aiai) and (b) (0-1(72) for ipQ = (strong coupling and 
anisotropic limit), O.Itt, 0.27r, 0.257r (no defect) , 0.37r,0.47r and O.Stt (free 
boundary condition). They are shown as a function of the distance r, in a 
log- log plot. 
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the sign of all spins on one side of the defect to change the sign of the defect 
strength. Thus our formulae again give the correct description of the system. 

From the agreement with the known cases, we believe our exact expres- 
sions ( ^.19|) and ( |5.2CI|) of the boundary two-spin correlation are correct, 



though we made some assumptions during the calculation. While the pa- 
rameter u is implicitly defined by the eq. (|5.7| ), there is no essential difficulty 
in solving the equation numerically. We used Mathematica to solve the equa- 
tion and then draw graphs of our results. We choose the location of the spins 
as shown in Fig. ^. We show the correlation functions for various strength of 
the defect in Fig. ^ as a function of the distance. In the short-distance limit, 
((TiCTi) converges to a single power-law which is independent of the defect 
strength. This is the expected bulk limit. In the large-distance limit, the 
correlation function is governed by another exponent, which depends on the 
defect strength. This is the boundary limit. In the next subsection, we will 
find the boundary exponent actually agrees with the boundary operator with 
the smallest dimension found in eq. ( |5.3|) . Our result interpolates between 
these two limits. 

In the "short-distance" limit, {(yiO'2) converges to a constant which de- 
pends on the defect strength. This actually corresponds to two spins located 
symmetrically about the defect line, and is related to the one-point function 
of the bosonic operator cosy?. In general, (o"i(T2) is smaller than (aicri) for 
the same defect strength and the same (horizontal) distance, as expected. 
Nevertheless, they asymptotically converges to the same power-law function 
in the large-distance limit for 6 7^ 0. Not only the exponent but also the 
constant prefactor is the same. This somewhat unexpected phenomenon will 
be further discussed in the next subsection using the boundary OPE. 

5.3 Boundary OPE of spin operators 



As discussed in Ref. ||T9[, once we obtain the boundary correlation function 
from the bulk OPE, we can calculate the boundary OPE from the boundary 
correlation function. 

The boundary OPE is the expansion about the boundary limit, where 
two operators are close to the boundary (real axis). Then the cross-ratio 
X, defined in ( |5.6| ), approaches 1. Thus we should expand the correlation 
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function in x = 1 — x. If we define u = e 



by 




(5.27) 



u is related to u = e 



by the modular transformation 



1 



(5.28) 



r 



The correlation functions (|5.19|) and (|5.20|) can be expressed in terms of 
u by the modular transformation of the elliptic theta function as 



From these expressions, we can read off the boundary OPE of ai and a2- 
In particular, for both (criO"i) and (o"icr2); the dimension of the intermediate 
boundary primary operators are given by 



where n is an arbitrary integer. They are exactly the boundary scaling di- 
mensions of the primaries in the "moving sector" (|5.2| ). Thus the "moving 
sector" is identified as the set of boundary operators generated by the bound- 
ary OPE of the spin operators. 

We note that, while the intermediate states in the bulk OPE are different 
for cTiO"! and o'ia2, the intermediate states in the boundary OPE are common 
between them. Only the OPE coefficients are different. Moreover, the OPE 
coefficient for the boundary operator with the smallest dimension is exactly 
the same between (cxicri) and (cti(T2). This means that the asymptotic behav- 
ior of these two correlation functions at large distance is exactly the same. 
(For negative b, they differ only by their sign.) The only exception is the 
free boundary case (6 = 0). In this case, there are two boundary operators 
for each boundary dimension and the two coefficients cancel exactly to make 
((Jicr2) identically zero. 




Afe = 2n^ + 



71 




(5.31) 
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The boundary OPE coefficients for ai±a2 are given by the sum of (or the 
difference between) those for ai and (T2. The smallest boundary dimension 
is given by 2(^)^ and 2(1 — respectively for cxi + (J2 and di — cr2- As 
discussed in Section pTT] , they are the only relevant boundary operators for 
this boundary condition. 



5.4 Correlation functions of the bosonic operators 

Here we consider correlation functions of the bosonic operators cos k(f. For 
k = 1, a/2 cos (y? is identified with the composite spin operator aia2 (at the 
same point). Thus the boundary two-point function of cos ip corresponds to a 
special case of the four-spin correlation function in the presence of the defect 
line: correlation function of two pairs, each of which consists of two spins 
located symmetrically about the defect line. 

The boundary one-point function can be determined by the method of 
Ref. PI as 

(cosM = ^^^- (5.32) 

In order to calculate the boundary two-point function {cos k(p cos kip) , we 
use a kind of mirror image method. Since it is symmetric under ip —>■ —ip, 
it is sufficient to consider only the fixed boundary condition ip = ipo at the 
boundary, apart from an ambiguity in the overall constant. The overall 
constant can not be determined by the mirror image method and will be 
fixed later by requiring the correct bulk limit. 

Given the boundary condition ip = ipQ at the boundary (real axis), the 
non-chiral boson field ip{z, z) can be represented by a single chiral boson (f){z) 
as 

ip{z,z) = ipQ + (t){z) - (t){z'), (5.33) 

where z' is the complex conjugate of z. Thus the boundary two-point function 
is mapped to the chiral four-point function as 

{coskip{zi,Zi) coskip{z2,Z2)) = 



1 ^^g±jfc[</3o+0(zi)-</.(4)]g±ifc|ipo+0(^2)-0(4)J^^ (5.34) 



^±ik[(po+(j>{zi)-(f){z[)]^±ik[ipo+<t){z2)-<l){z!2)] 

where it's are summed independently. 
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Using the standard result on the correlation functions of the vertex oper- 
ators 

^ -Q^^^. _ ^^^.,a,/4 (J^a, = 0), (5.35) 

j j<k j 

this reduces to 

C 2 2 

{cosk(p{zi,zi) coskip{z2,Z2)) = T- VT277[a;~'' + x'' ^^cos2kipo] (5.36) 

where C is a constant and yi, ?/2 and x are defined similarly with eq. (|5.19|) . 
Requiring the correct bulk limit 

(cos ki{>{zx, %) cos ki{>{z2, ^2)) 771 1.2/4 , (5-37) 

l\Z\ — Z2\ ' 

C is fixed to be 1/2. 

Let us compare our result ( p.36| ) for = 1, which is the special case of 
four-spin correlation function, with known results. For periodic {ipQ = vr/4) 
and antiperiodic {ipo = 37r/4) boundary conditions, our result agrees with the 
bulk four-point spin correlation of the Ising model for the special location of 
the spins. For the free (y^o = ^/2) boundary condition, our result agrees with 
the square of the boundary two-spin correlation for the free boundary 
condition. Again we found consistency of our approach with the known 
results. 



5.5 Two-spin correlation near the end of the defect 
Une 

So far we have discussed only the defect line without ends (infinitely long or 
periodic). Let us consider the neighborhood of the end of the defect line. We 
discuss only the correlation function of two spins located symmetrically about 
the defect line, as shown in Fig. ^. This can be reduced to the one-point 
function (a/2 cosy?) = (o"icr2) after the folding, with a changing boundary 
condition. 

We regard the end point as the origin, and make the conformal mapping 
z = e"^ to the strip with different boundary conditions at two sides. Then 
the one-point function is given by the vacuum expectation value with two 
specified boundary conditions. The Hilbert space of the strip is classified 
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spin operator 



defect line 




end of 
the defect 



Figure 5: The correlation function of two spins located symmetrically about 
the defect line. It is characterized by the distance from the line y and the 
angle 9. 



47 



into several sectors by the zero-mode. We only need to consider the lowest 
energy zero mode to calculate the vacuum expectation value. Thus we fix 
ip = ipo at the boundary corresponding to the defect and ip = 7r/4 at the 
other boundary. 

We employ the mirror image method used in the last subsection. In order 
to ensure the boundary conditions on the two boundaries, we first need two 
mirror images as 

p>{w,w) ~ ^ + (^0 - 7)- + 4>{w) - 4>{w') - 4>{w' + 2m), (5.38) 

where is a chiral boson, w' is the complex conjugate of w and v = Imw. 
However, the mirror image about one boundary confiicts with the boundary 
condition on the other boundary. Thus we need the mirror image of the 
images. In this way, we need infinitely many mirror images as 

(f{w, w) = — + {(fo ~ ^) X! + 2vrm) + ^ (j){w + 27rm). (5.39) 



Using eq. (|5.35| ), the one-point function is given by 



/ 00 n \ 1/4 



(cosv?) = const, cos [- + (v^o - 7)-] H 



4 A'tt' \J':^2v + 2Tm^ 



cos [f + {ifo 



const. V,_xi/4 • (5.40) 



smf 



Going back to the half plane, we obtain the one-point function, or the special 
two-spin correlation before folding as in Fig. ^, near the end of the defect 
line. It is given by 

(.i(...>a(---°)>-V2 °°'''^^^^°;'''' , (5.41) 

where z = is the end of the defect line, y = Imz, 9 = arg z, and the overall 
normalization is determined from the ^ — >^ 0, tt limit. In the ^ — > limit, the 
spins are located far apart from the defect line and the correlation function 
should approach the bulk value {2y)~^^'^. In the 6 ^ n limit, the spins are 
located symmetrically about the defect line and are far apart from the end 
of the defect. Thus the correlation function should converge to that for the 
infinitely long defect ( p.32| ). 
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5.6 Free energy of the finite length defect line 



Here we discuss the free energy due to the presence of the defect hne of 
the continuous Dirichlet universahty class ( p.36|) . A defect hne with a finite 
length corresponds to a finite interval with a different boundary condition, 
after the folding. The partition function in the presence of the defect line 
is expressed by the two-point function of a boundary condition changing 
operator [^. (For further discussion, see Appendix.) 



The dimension of the boundary condition changing operator is determined 
by the ground-state energy of a strip with boundary conditions corresponding 
to the defect line and no defect (6 = 1) on the two sides. From the partition 
function ( p. 33 ), the dimension is found as 



Ab = - 
2 



vr 4 



(5.42) 



The universal term in the free energy of a defect of length / on an infinite 
plane is thus given by 



i^univ(0 = -kT log Z = kT Ablogl. 



(5.43) 



The universal free energy of the defect line along the circumference L of a 
cylinder can be obtained by a conformal mapping z = tan tcw /L as 



F,„i,(/,L) = fcTAfelog 



L Til 

— sm — 

77 L 



(5.44) 



However, the above is not the only term in the free energy. There are 
also non-universal terms. Firstly, there is a non-universal shift of the ground 
state energy depending on the boundary condition, as seen in the second 
term in eqs. ( 2.20 ) and ( 2.21 ). This is due to the phase shift of the fermion 
momenta 



, which depends on the fermion energy. This effect, which is 
known as Fumi's theorem ||3^ in the general case, depends on the dispersion 
relation and the phase shift in the entire Brillouin zone, and thus is non- 
universal. In general, we expect a shift of the ground-state energy of order 
1/a and the free energy proportional to the defect length as Al, where A is a 
constant of order 1/a. An exception is the antiferromagnetic defect (h = —1). 
In this case, the quantization of the wavevectors is exactly the same as in 
the no defect case and only the parity assignment is reversed. Therefore the 
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non- universal term proportional to the length is absent {A = 0). This is also 
required from the fact that the antiferromagnetic defect is exactly the two- 
point correlation of disorder operators , which decays purely algebraically. 

Secondly, the free energy of the defect line should vanish when there is 
no defect (/ = 0). Hence a constant term of order logo, where a is the 
cutoff length (lattice constant) is required. The combined result would be 
approximately kTA^ log -, which is only valid for I » a. The constant term 
of order log a is non-universal and can be regarded as an effect of the ends 
of the defect line (or, of the insertion of the boundary condition changing 
operator.) Moreover, the dependence of the free energy on I when I ~ a 
should be non-universal. If the defect length approaches to the circumference 
of the cylinder {I L), the dependence on / becomes again non- universal for 
[L — 1) a. When / = L, however, the defect line completely wraps around 
the cylinder; there is no longer a non-universal constant due to the ends of 
the defect line. Nevertheless, a universal constant of order 1 still remains 
in general. After the folding, we have a new single boundary condition. 
The universal remaining constant is nothing but the contribution from the 
"ground-state degeneracy" ||3^: — kT log g. In the present problem, (7 = 1 for 



any continuous Dirichlet boundary state, as will be discussed in Section |^ 
Thus the universal constant for / = L is zero. 

As discussed above, the non-universal part in the free energy of the defect 
line is rather large except for h = —1; it may be difficult to measure the uni- 
versal part in (numerical) experiments. (The dimensions of boundary condi- 
tion changing operators are measurable in quantum impurity problems |[40|| . ) 

5.7 Duality and disorder correlation 



Our model ( |2.1| ) has a defect line where the strength of vertical (i.e. or- 
thogonal to the defect line) links are altered. Under the duality transfor- 
mation our model ( |2.1D transforms to the Ising model with the 



defect line where the strength of the horizontal coupling is altered. This 
"horizontal" model was studied by several authors 0, 3. The duality 
transformation exchanges the order and disorder operator. It is illuminat- 
ing to discuss the duality transformation in the quantum model under the 
r-continuum limit. Our "vertical" model corresponds to the quantum Hamil- 
tonian (|2.16| ). The duality transformation is introduced by defining the 
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new operators /t as in 



fi [n) 
ff{n) 



cr^(n)dr^(n + 1), 



(5.45) 



where a's are defined as in ( |2.16| ). ft^ is the disorder operator. 

We apply this dual transformation for — oo < n < oo. Using the new 
operators ft, the Hamiltonian (|2.1(j|) is expressed as 



H 



n 



E 



(5.46) 



Now the Hamiltonian is again the transverse Ising model, but the only the 
strength of the transverse field is altered at n = 0. This corresponds to the 
variation of the horizontal link strength. 

Let us consider the model by a perturbation from the model without the 
defect {b = 1) as in Ref. 0. In the continuum limit, the variation of ei- 
ther vertical or horizontal coupling correspond to the perturbation by the 
same energy density operator. Thus, if the continuum perturbation theory 
is valid in principle, the "horizontal" model is equivalent to the "vertical" 
model (|2.1|) under some change of the parameter y^o- Assuming this, the 
Dirichlet boundary state is transformed to the Dirichlet boundary state with 
a different parameter ip^ under the duality transformation. The transforma- 
tion rule for the parameter ipo can be determined from Brown's relation 
between the surface exponents of the order /disorder fields: 



2A,(a) + V2Ab(/x) = 1. 



(5.47) 



This relation together with our result ( p. 3D gives the transformation of ipo 
under the duality transformation: 



^0 ^0 



f-^o (0<<^o<f) 
<^o - f (f < V'o < vr) 



(5.48) 



Thus the correlation functions of two disorder operators in the Dirichlet 
boundary state are given by 



(5.49) 
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6 Properties of the continuous von Neumann 
universality class 

Here we discuss the correlation functions near the defect hne that is identified 
with the continuous von Neumann boundary state ( 4.34|) . We also identify 



the corresponding defect line in the one-dimensional quantum model. 
6.1 Two-spin correlation functions 

We can calculate the two-spin correlation functions {(JiCTj) in a similar manner 



to what we did in Section |5]^. The only difference is the boundary state, and 
consequently the coefficients defined in (|5.9| ). The boundary state ( [4.34| ) 
determines the coefficients as 

AnV2,A) ^ f 1 = 0) . . 

^ \ V2cos(2n(^o) (n> 1) ^ ' 

Firstly, it is easy to see that (criO"2) vanishes for any value of ipo, since the 
bulk OPE (Ti and cr2 only generates the operators as in ( ^.141) . 
The bulk OPE coefficients of ajaj (j = 1, 2) are given by 

a,.,(nV2,A) = ^[Y^) ' (6-2) 

/ 1 \"V2 

C^,.,iny2,A) = (-1)"V2(t7j) • (6-3) 



16 



Thus the correlation functions are given by 



{ai{zi,zi)ai{z2,Z2)) 



1/ 



1 CO „ ^ 

t -sr^ .(nV2,A)^ 



2 



J2^B ' C'ai<7i(n2/2,A)/(l7, T^, 



1/8 



[1 + ^^2 cos {2n(po)u 4 



Me''^\V^), (6.4) 



1 V^' 1 



{a2{Zi,Zi)a2{z2,Z2)) 
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\4:yiy2xJ i^sl^j^o 2 16 

/ 1 \^^^ 1 °° 2 

\4:yiy2xJ Mu) 

/ 1 \'/' 1 

= I ^r^Me''^',V^). (6.5) 

They are no longer symmetric under the exchange of ai and (T2, or equiva- 
lently the reflection about the defect line before the folding. We notice that 
these correlation functions are actually identical to those for the continu- 
ous Dirichlet universahty class (|5 .191 , 15.201 ) for an appropriate value of the 
parameter. We summarize this observation below 

(criO-l)jV(<^o) = {(^1(^1) D{vo=vo) 

{cr2cr2)N{^o) = (criO-i)z)(¥>o=V2-w) (6-6) 

(o"lO"2)Af(<^o) = 0) 

where () N{>f>o) and () d((/3o) correlation function in the von Neumann and 

Dirichlet boundary states, respectively. It is also possible to calculate other 
quantities discussed in Section |^ for the continuous von Neumann universality 
class. 



6.2 Identification of the defect line 

Naively the result {o'ia2) = implies that the two layers are completely 
decoupled. If they are completely decoupled, the boundary state must be 
Iff) or equivalently the Dirichlet one with = tt/2. However, for example, 
at <^o = 7r/4 the correlation functions (cticti) and {(72(^2) behaves exactly as if 
there is no defect, while (cri(T2) still vanishes. These behaviours are definitely 
different from those in \ f f) boundary states. Hence there should be some 
sort of coupling between two layers. 

We have not yet found an appropriate defect line in the classical model 
which corresponds to the continuous von Neumann boundary state. However, 
we have succeeded in finding an appropriate point defect in the corresponding 
one-dimensional quantum model. It is given by the Hamiltonian: 

H = -J2 ^"H - E ^'(^ - - &'7'(-1)<t"(0), (6.7) 
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where the operators are defined similarly as in eq. ( p.l6| ). 



We show that this model is related to the previous model ( 2.16 ) by a 
duality transformation ( p. 451 ) on one side of the defect only. The duality 
transformation on the half-line n > is defined by 



0<m<n 

= ci\n)a\n + l). 



(6.8) 
(6.9) 



The new operators /t's are defined only for n < and commute with a's at 
n < 0. They obey the same commutation relation as the original a's do. The 
Hamiltonian (16. 71) is rewritten as 



H 



- E - l)a\n) - ^ /i^(n - l)fi^{n) 

7i<0 n>0 

-ha'{-l)fi'{{]). 



(6.10) 



This is exactly the Hamiltonian ( 2.16|) , under the redefinition fi —>■ a. 

In this way, the correlation functions in our model {)m are identified with 
the correlation functions in the Dirichlet boundary condition as follows: 



{(^icri)M 



(CX10"2) 



M 



(o-iO-i)d(<^o), 

(Ai2/i2)D(</30) 
((Xl/i2)z)(vPo) = 



(6.11) 

(6.12) 
(6.13) 
(6.14) 



where n is the disorder operator ||S^, (p^ is defined in ( |5.48| ) and (po is deter 



mined by the parameter b as in (|3.37| ). 

These results exactly agree with the correlation functions ( |6.6| ) in the von 
Neumann boundary condition if ipo = (pQ. Therefore our model (|6.7| ) pre- 
sumably corresponds to the continuous von Neumann boundary state ( [4.34| ) 
in the continuum limit. 
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7 Renormalization group flow among the bound- 
ary states 

7.1 The relative stability of the boundary states 

Now we discuss the relative stability and the renormalization group flow 
among the boundary states. In this section, the orbifold radius r is again 
fixed to 1 as we discuss the Ising defect problem. The universal "ground-state 



degeneracy" g , which is given by the inner product of the boundary state 
and the ground state, turns out to be 1 in the continuous Dirichlet boundary 
state ( [4.32| ) for any value of (^q. Similarly, g = ^J2 in the continuous von 
Neumann boundary state ( [4 .341 ) independent of (p^. The number g plays 
the role similar to the central charge c in the bulk conformal field theory. 



According to the "gf-theorem" |^ , g decreases along a renormalization group 



trajectory connecting two conformally invariant boundary conditions. 

In ( |4.32|) and ( [4.34| ), the boundary state depends continuously on the de- 



fect strength 6, namely there is a fixed line of conformally invariant boundary 
conditions. Thus it is natural that g is constant along two continuous line of 
boundary universality class, as central charge c is constant along a fixed line 
of conformally invariant bulk theories. Namely, no renormalization group 
fiow can occur along these lines. The values of g in the 8 discrete boundary 
states are simply the product of g for the two Ising layers. It is known 
that (? = 1 for the free Ising boundary condition and g = for the fixed 

Ising boundary condition. Thus g is 1/ ^J2 for the endpoint Neumann bound- 
ary states (I t /) etc.) and 1/2 for the endpoint Dirichlet boundary states 

(ITT))- 

By the "^f-theorem" , the continuous Neumann boundary state is the most 
unstable among these boundary states. The next most unstable one is the 
continous Dirichlet, then the endpoint Neumann, and the most stable ones 
are the endpoint Dirichlet. If there is no restriction on the defect line, it is 
thus generically attracted to one of the endpoint Dirichlet boundary states. 
In order to realize the other universality classes, we must protect the defect 
line from the relevant operators, for example by an appropriate symmetry. 
From the expression of the boundary states as the tensor products of the 
Ising boundary states, it is clear that the discrete endpoint boundary states 
break the global symmetry: o —a. Thus if we do not break the global 
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^2 symmetry, the defect line is renormalized to one of the continuous states. 
On the other hand, the endpoint Neumann boundary states are expected to 
have a different Z2 symmetry: Ising spin flip only on the one side of the 
defect line. For example, | t /) is invariant under the Ising spin flip in the 
second layer. Therefore the defect line will be renormalized into one of the 
endpoint Neumann ones if we break the global Z2 symmetry but preserve 
the other symmetry. 

Analysis of the operator contents is useful to further investigation of this 
problem. We can obtain the spectrum of the boundary operators in the 
continuous Neumann boundary states in a similar manner as in Section 
The partition function on the strip with the same Neumann boundary on 
both sides is given by 

-i 00 -1 00 

where q is defined as in ( |2.26| ). 

From this expression, we can read off the whole spectrum of boundary 
operators. In particular, the dimensions of the relevant operators (other than 
the identity) are given by 1/2 (twice), 2(^)2, |(1 - '^f. From the two- 
spin correlation functions ( |6.6[ ), the latter two are identified with the spin 
operators ai and a2- The remaining two with dimension 1/2 are identified 
with 0"icr2- The multiplicity 2 is understood as follows. Since the dimension 
1/2 is independent of the parameter (pQ, we may consider the simplest case 
6 = 1 (<^o = vr/4). At this point, the problem is reduced to the Ising model 
without any defect after the duality transformation. Now we can unfold the 
system and reduce the boundary OPE between ai and (72 to the bulk OPE 
between a and /i, which is given [^, as 



a{w, w7)/i(0, 0) ~ w^'^w-'^^^i) + w-^/^u7^/V- (7.2) 

Thus there are actually two operators if) and generated by the OPE of a and 
/i. The coefficients depend on the relative location of a and /z. Translating 
to the folded picture, two boundary operators are generated by the OPE 
of o"i and o"2 and the relative coefficients depend on the location of ui and 
(72. While the above argument is only valid for the special point h = 1, it 
is natural to expect that two boundary operators are generated by the OPE 
for any value of h (or (po). 
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Now let us turn to the stability of the Neumann boundary state. Since the 
spin operators ai and a2 are relevant, the perturbation by these operators will 
drive the system away from the von Neumann boundary state. The resulting 
state will be one of the discrete boundary states. For example, if we perturb 
by cTi, the boundary will be presumably renormalized to | f /) (or | i /)) 
boundary state. The remaining relevant operators are generated by aiO'2- 
Under this perturbation, the boundary will be presumably renormalized to 
the Dirichlet boundary state with some value of <^o- These conclusions are 
consistent with our discussion based on the "(yf-theorem" . 

To achieve the continuous Neuman boundary state in the infrared limit, 
we must protect the boundary from these relevant operators. The model (|6.7| ) 
has the following Z2 x Z2 symmetry: 

a^{n) -a%n) {n < 0) with a^(0) -(7^(0) , . 

a'ln)^-a'ln) (n > 0) ^ 

In the continuum limit, it may be simply interpreted as 

(Ji — > —0"! and (T2 —(72. (7.4) 

Any of the relevant operators ai, (J2 and o"icr2 breaks this Z2 x Z2 symmetry. 
This means that the Z2 x Z2 symmetry protects the von Neumann boundary 
state. However, the Dirichlet boundary state with ipo = 7i/2 (i.e. |//)) also 
has the Z2 x Z2 symmetry. Hence a defect line with the Z2 x Z2 symmetry may 
not renormalize into the von Neumann boundary state; an obvious example 
is our original model ( p.l|) with 6 = 0. 

As we have discussed, the aia2 perturbation will presumably drive the 
defect line to the Dirichlet boundary state, where only the single global Z2 
symmetry ai^2 — cri,2 exists. This is natural since cti(T2 breaks the Z2 x Z2 
symmetry to Z2. On the other hand, the reflection symmetry ai ^ o"2; which 
was absent in the continuous Neumann boundary state, is recovered in the 
continuous Dirichlet boundary state. 

The boundary operator content in the continuous Dirichlet boundary 
state is given in Section The two boundary operators with the dimensions 
2{(Pq/tt)^ and 2[{n — (Pq) / n]"^ are only the relevant boundary operators. As we 
have shown in Section they correspond to ai±a2- At the "periodic" point 
(fio = 7r/4, Af) = 1/8, 9/8 which correspond to a and da/dx ~ (ai —a2)/a. At 
the "free" point ipo = n/2, At, = 1/2, 1/2 since cxi and a2 are independent. 
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Thus, if we preserve the global Z2 symmetry 0\^02 — — o"2, the fixed 
line is stable. If we perturb the defect line by the relevant operator ui ± (T2, 
the boundary state will be presumably renormalized to one of the magnetized 
boundary states (| It) ^tc.) Such a fiow is established at the "free" point 
{(fQ = vr/2), where the boundary condition of the two independent Ising 
models are renormalized. We expect the renormalization group fiow from 
the entire fixed line to the discrete magnetized boundary states, consistently 
with our arguments based on "^f-theorem" . 



8 Generalized defects and universality 



Grimm considered generalizations of the defect line. Here we discuss his 
result from our general consideration of defect universality classes. 

His first ( "integrable" ) generalization of the defect line is given by the 
quantum Hamiltonian 



N-l 

E 



(8.1) 



where 



= cos (ha + sin (ha 



y 



(8.2) 



(Eq. ( ^.11) corresponds to eq. (3.1) in Ref. but we have changed the 
notation.) When = 0, this reduces to the "ordinary" defect described by 
eq. (|2.16| ). Grimm obtained the exact spectrum of this model. Apart from 
the non-universal shift in the ground-state energy, the result was identical 
to the "ordinary defect" , i.e. continous Dirichlet boundary state ( [4.32| ) with 
some value of ipo. He noticed that the model (|8.1| ) preserves the global Z2 
symmetry (simultaneous rotation about a;-axis by angle vr at each site) and 
thus it belongs to the same universality class as the "ordinary" one (|2.16|) . 
These are completely consistent with our arguments based on "^f-theorem" 
and the boundary operator contents. 

He further introduced another generalization given by 



N-l 



1 



[a^ + Pa^^{cosipal - sin-^cri^) 



(8.3) 
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where < ip < it. (This corresponds to eq. (4.1) in Ref. |T^ but again we 
changed the notation.) This model is not exactly solvable and he obtained 
the spectra numerically. He calculated the difference of the universal term in 
the ground-state energy between (|8.3| ) with P > and eq. (|2.16| ) with 6 = 0. 
The numerical result was consistent with his conjecture that the difference 
is given by 0,1/16 and 1/2 (in unit of 27rv/N), respectively for < 7r/2, 
ip = 71/2 and ip > n/2 (independent of /3 as long as /5 > 0.) He argued that 
this means that there are three universality classes for this model. 

Let us reconsider his model (|8.3|) from our viewpoint. As he noticed, 
this model no longer keeps the global Z2 symmetry. From our arguments 
in the previous subsection, we expect that the defect line will generically 
renormalize into one of the endpoint Dirichlet boundary states. Actually, we 
can understand as follows. (Here we restrict to the case in which (3 > 0.) 
We expect (cr^) > due to the transverse field in the Hamiltonian ( |8.ij| ). 
Thus (o"^) is expected to be negative for all < ip < n, due to the effective 
field generated by af. Then af is affected by the effective field generated 
by cr^. Depending on the sign of the coupling, (af) should be negative for 
< Ip < n/2 and positive for7r/2<?/)<7r. At ip = n/2, neither direction is 
favored for af. 

Thus we expect that the corresponding boundary state is | ||), | | /) 
and I It), respectively for < ■i/' < n/2, ip = n/2 and n/2 < ip < n. From 
the discussion of the previous subsection, the endpoint Neumann boundary 
state I I /) is unstable and would be realized only at = n/2. Actually, at 
this point we have the extra symmetry: —a^ only on one side of the 

defect. This symmetry protects the endpoint Neumann boundary state, as 
we argued. Changing ip from n/2 corresponds to the relevant perturbation 
and the defect line will renormalize to the endpoint Dirichlet: | \\) or | ||) 
depending on the direction of the perturbation. We also note that, although 
his model ( ^.3|) oXtp = n /2 looks similar to our model (|6.7|) which corresponds 
to the continuous von Neumann boundary state, they are actually different. 
The transverse field at site ra = exists in the former, breaking the x Z2 
symmetry ( [7. 3] ) to the single Z2 symmetry. 

Grimm calculated the ground-state energy of the model (|8.3| ) on a circle. 
In our formalism, it corresponds to a strip of width A^/2, with the defect 
boundary state on one side and the "periodic" (Dirichlet with (/jq = 7r/4) 
boundary state on the other side. The partition function of the system is 
readily derived; only the untwisted sector of the defect boundary states con- 
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tributes to the partition function. From the partition function, the ground- 
state energy of the strip is given by 1/16, 1/8 and 9/16 in the unit of 2nv/N, 
respectively for | ||), | | /) and | |t). Following his definition, we subtract 
the ground-state energy in the case the defect boundary state is |//), to 
obtain 0, 1/16 and 1/2 respectively for those three cases. Thus our analysis 
is in complete agreement with Grimm's. Furthermore, his conjecture based 
on c = 1/2 character of the Virasoro algebra also agrees with our boundary 
states which is tensor products of Ising boundary states. 

The "new" universality class of the defect lines he obtained was actually 
equivalent to boundary conditions of two independent Ising models. Never- 
theless, his general discussion on the universality of the defect lines seems 
quite correct. 

Recently, Karevski and Henkel discussed another kind of generaliza- 
tion. They studied a junction of S" = 1/2 and S = 1 quantum Ising chains. 
In the continuum limit, the bulk universality class of the S = 1 Ising chain 
is identical to that of standard = 1/2 Ising model [|^. Thus, when both 
chains are critical, we expect that the universality class of the junction is 
identical to that of the defect discussed in this paper. Actually, they numer- 
ically studied [44| the junction with the Ising Z2 symmetry, and found the 
universal behaviour corresponding to the continuous Dirichlet boundary con- 
dition in our terminology. This also confirms that the notion of universality 
applies to the defect problem. 



9 Summary and Discussion 

We studied the continuum limit of the two-dimensional critical Ising model 
with a defect line. Folding the system at the defect line, we map the system 
to the critical Ashkin- Teller model [^, 3^ at the decoupling point with 
boundaries. Thus we can classify the universality classes of the defect lines 
by the boundary states of a conformal field theory. 

Based on the exact partition function, we identified the boundary state 
corresponding to arbitrary strength of the defect, in terms of the Ashkin- 
Teller Ishibashi states. Furthermore, we discussed the possible boundary 
states in the Z2 orbifold of free boson, which is identified with the continuum 
limit of the critical Ashkin- Teller model. There are two continuous lines of 
boundary states, and eight discrete ones for general values of the compact- 
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ification radius r. The special case r = 1 is applied to our problem of the 
defect lines in the Ising model. It is shown that the obtained orbifold bound- 
ary states include all the known defect lines, and one new one-parameter 
family universality class of defect lines corresponding to the orbifold von 
Neumann boundary state. We conjectured that the above boundary states 
exhausts all the possible boundary states on orbifold for generic value of 
r; if this is true, we have identified the complete set of universality classes of 
defect lines in the Ising model. 

From the boundary state, we obtained the complete spectrum of boundary 
operators. We calculated the exact boundary two-spin correlation functions 
for arbitrary strength of the defect, employing the method of Cardy and 
Lewellen and Zamolodchikov's solution [^] for c = 1 conformal block 
of spin operators. We also calculated the boundary two-point functions of 
bosonic operators, one of which is identified as a special case of the four-spin 
correlation function. In addition, we obtained a special case of the two-spin 
correlation function near the end of the defect line, the universal term in 
the free energy of the defect line and the correlation function of the disorder 
operators. While we made some assumptions during the calculation, our 
approach, which satisfies several nontrivial consistency check, appears to be 
valid. 

We also calculated the two-point spin correlation functions in the new 
one-parameter family of defect lines, and identified the corresponding quan- 
tum model using the duality transformation. For the 8 discrete universality 
classes, calculation of the correlation functions reduces to that in the Ising 
model with a boundary, which is already solved [0. Thus we have obtained 
two-point spin correlation functions for all defect universality classes found so 
far (i.e. for all the universality classes, if our conjecture is correct.) Finally, 
we discussed the renormalization-group flow among the defect universality 
classes based on the "(/-theorem" and the operator contents. We explained 
Grimm's results on the generalized defects, applying our observations. 

A possible extention of our present work is the analysis of the boundary 
states in the general critical Ashkin- Teller model, i.e. for general value of the 
compactification radius of the orbifold. While two continuous families and 
eight discrete boundary states exist in general, their relative stability will 
vary. Moreover, for special values of the radius, additional boundary states 
may appear, as in the case of the free boson ||43[| . 

In Ref . H] , the first attempt was made to apply boundary conformal field 
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theory to the defect problem via folding. However, it was essentially repro- 
duction of the earlier results. On the other hand, our present result provides 
a more detailed description of defect lines in the Ising model, than that in the 
earlier literatures. This demonstrates the power of boundary conformal field 
theory for the application to the defect problems. In particular, the solution 
indeed posesses c = 1 structure which is absent in the original c = 1/2 prob- 
lem before the folding. Thus the folding seems a really necessary procedure 
to treat the defect problem in terms of conformal field theory. From the gen- 
eral viewpoint on boundary conformal field theory, the present model seems 
one of the simplest nontrivial models with an infinite number of boundary 
conformal towers. Cardy's original treatment p|, ^ mainly focused on min- 
imal models which have a finite number of conformal towers. Most of the 
successful applications to c > 1 have also been to theories with a finite 
number of primaries with respect to current algebra. When the number of 
primaries cannot be reduced to a finite one, there are several technical diffi- 
culties [|]. In spite of those difficulties, we found that the present problem 
is tractable, though we have not yet completely understood the structure 
of the theory as we do in minimal models. We hope that the observations 
on the present model will be useful in the further development of boundary 
conformal field theory with infinitely many conformal towers. 
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Appendix: free energy and the boundary con- 
dition changing operator 

By folding, a finite length defect line is mapped to a changed boundary 
condition for a finite interval. Let us discuss the free energy associated to 
such a finite interval of changed boundary condition. Using the transfer 
matrix (Hamiltonian) in the direction parallel to the boundary (defect line). 
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the partition function of the system may be written as 



Z{1) = Tre 



-{L-1)Ho^-IHd 



(A.l) 



where / is the interval of the changed boundary condition (length of the 
defect), L is the system size in the direction parallel to the boundary and Hq 
and Hd are Hamiltonian for the original and changed boundary conditions. 
The partition function in the absence of the defect line is 



Zo = Tre-^^o. 



(A.2) 



On the other hand, a change in the boundary condition can be associated 
with a "boundary condition changing operator" ||28|. The boundary con- 



dition changing operator acts in an enlarged Hilbert space with all possible 
boundary conditions, and maps a whole subspace with a particular boundary 
condition into each other. By the boundary condition changing operator U, 
the Hamiltonian is mapped as 



U^HM = H 



D 



(A.3) 



In their pioneering work, Schotte and Schotte introduced such an oper- 
ator for the free boson. (See also [0.) Here we discuss the non-orbifold free 
boson. The boundary condition changing operator that maps the Dirich- 
let boundary condition = to the Dirichlet boundary condition with a 
different boundary value <y9 = 5 is given as 



U 



exp 



(A.4) 



where <f{0) is the quantum operator corresponding to the (p at the boundary. 
Here the operators are quantized, regarding the direction parallel to the 
boundary as a imaginary time. As they claimed, U is apparently a unitary 
operator. Then it follows 



,-lHo 



u 



,-IHd 



and thus 



Z{1) = Tr 



(A.5) 



(A.6) 
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Namely, the ratio of Z{1) and Zq is identified with the two-point correlation 
function of the boundary condition changing operator: 



^ = {u^mm, (A.7) 

where the expectation value is defined in the original theory with the Hamil- 
tonian Hq. In the following, we assume the thermodynamic limit L —>■ oo. 
The asymptotic behavior for large / is governed by the lowest dimension 
of the boundary condition changing operator. This gives 

{U^{l)U{0)) ~ const./-^" (A.8) 

The free energy of the changed boundary condition (defect line) of a finite 
length is given by taking the logarithm as 

Z(l) 

F{1) = -kTlog-^ = const. + kT Ablogl. (A.9) 

The constant term is non- universal. From the dimensional analysis, we ex- 
pect it to depend on the short- distance cutoff (lattice spacing) a as roughly 
— log a In general, there can be also a non-universal shift in the energy for 
a different boundary condition. This contributes extra Al, where A is a 
non- universal constant of order 1/a, to the free energy. 

However, there is a subtlety in the above argument: U cannot be a unitary 
operator in the standard sense, since Ho and Hd have different spectrum in 
general. This peculiarity is presumably related to the presence of infinite 
degrees of freedom in the field theory. Thus the identity ([A. 51) and the above 



argument may be questionable. We argue that U is still a unitary operator in 
some weak sense and the above discussion on free energy is valid in general. 
For the free boson case, we can demonstrate the validity by calculating both 
sides of the equation (|A.7|) independently. 



Let us consider the free boson with the Lagrangian ( |3.2| ). We take the 
boundary as x-axis {y = 0) and assume that the system lies in the upper 
half plane. We consider the Dirichlet boundary condition = 6 for the 
interval < x < I, and = elsewhere. The corresponding boundary 
condition changing operator is given in (|A.4|) . The asymptotic form of the 
two-point correlation function of the boundary condition changing operator 
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is given by ( |A.8| ) with A;, = ^^/vr^. On the other hand, the Free energy in 
the absence of the changed boundary condition is given by 

Zo = J P(/pe-5^/'^'"(^'^^)', (A. 10) 

where (p = aX the boundary {y = 0). In the presence of the changed 
boundary condition, we write 



(f = iPc + iPg, (A. 11) 



where ipc satisfies 



d'cpc = 

/ ^\ f (X < 0, / < X) / A 1 r^N 

^c(x,0) = 1^ [q^;^^^ (A.12) 
fcix,y) (y->oo) 

(fg is only required to satisfy = at the boundary, which is exactly the 
same condition ip obeyed in the absence of the changed boundary condition. 
In other words, the change in the boundary condition is completely absorbed 
by ipc- The partition function in the presence of the changed boundary 
condition is given as 

Z{1) = J P(^ge-^/'^'"[^^(^«+^^)l'. (A.13) 

The action can be decoupled as 

1 



2tt 



d'xid.if.y + id.ipgY, (A.14) 



since the cross term is shown to vanish by partial integration. Actually there 
is an infinite number of solution for ip^ and we must sum over such allowed 
ifc, apart from the integration over (pg. However, in the large-/ limit, only 
the one with minimum action dominates. For \S\ < vrr, ipc with minimum 
action is given by 

(Pc = - ftan"^ - tan"^ -) . (A.15) 

TT V X — I Xj 
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(The ones with larger value of the action is given by replacing 6 hj 6 + 2iTr, 
etc.) The leading term in the partition function is given by 



Z{1) = J Vifgexp [-^ J d^x{df,!fcf + {df,fq 



The functional integral over ipq is identical to Zq in eq. ( A.IO ). The action 



of (pc can be calculated using the explicit form (|A.15|) . Actually it diverges 



at the singularities {x,y) = (0,0) and (/,0). Introducing the short-distance 
cutoff a, we obtain 

-L f£^{d,^^f = ^hg-. (A.17) 

Thus we finally get 

^ = {jf"", (A.18) 

which agrees with the two-point correlation function of the boundary condi- 
tion changing operator. Thus the validity of the relation ( [A.7| ) is explicitly 
shown in the free boson case. 

This calculation can be generalized to the orbifold boson. On the other 
hand, in general conformal field theories, it is difficult to calculate directly 
the partition function in the presence of the changed boundary condition. 
We expect that the relation ( |A.7| ) also holds in such general cases. 
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